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ABSTRACT 

Context. The [C n] 158 pm line is one of the dominant coolants of the ISM, and an important probe with which to study the star for¬ 
mation process. Recent Herschel/HIFI and SOFIA/GREAT observations showed that assuming the total velocity-integrated intensity 
of this line is directly associated with the star-forming material is inadequate. 

Aims. We probe the column densities and masses traced by the ionized and neutral atomic carbon with spectrally resolved maps, and 
compare them to the diffuse and dense molecular gas traced by [C i] and low-/ CO lines toward the star-forming region M17 SW. 
Methods. We mapped a 4.1 pc x 4.7 pc region in the [C i] 609 pm line using the APEX telescope, as well as the CO isotopologues 
with the IRAM 30m telescope. Because of the velocity-resolved spectra, we analyze the data based on velocity channel maps that 
are 1 kms -1 wide. We correlate their spatial distribution with that of the [C n] map obtained with SOFIA/GREAT. Optically thin 
approximations were used to estimate the column densities of [C i] and [C n] in each velocity channel. 

Results. The distribution of the emission from the isotopologues 13 CO, C 17 0, and C 18 0 resembles more closely that of the [C i] 
emission than that of the 12 CO emission. The spatial distribution of the [C i] and all CO isotopologues emission was found to be 
associated with that of [C n] in about 20%-80% of the mapped region, with the high correlation found in the central (15-23 kms -1 ) 
velocity channels. 

Conclusions. The excitation temperature of [C i] ranges between 40 K and 100 K in the inner molecular region of M17 SW. Excitation 
temperatures up to 200 K are found along the ridge. Column densities in 1 km s -1 channels between ~10 15 cm -2 and ~10 17 cm -2 were 
found for [C i]. Just ~20 % of the velocity range (~40 kms -1 ) that the [C n] line spans is associated with the star-forming material 
traced by [C i] and CO. The total (integrated over the 0-40 kms -1 velocity range) gas mass estimated from the [C n] emission gives 
a lower limit of ~4.4xl0 3 M 0 . A very large fraction of at least 64% of this mass is not associated with the star-forming material 
in M17 SW. We also found that about 36%, 17%, and 47% of the [C n] emission is associated with the H n, H I, and H 2 regimes, 
respectively. Comparisons with the H41cr line shows an ionization region mixed with the neutral and part of the molecular gas, in 
agreement with the clumped structure and dynamical processes at play in Ml7 SW. These results are also relevant to extra-galactic 
studies in which [C II] is often used as a tracer of star-forming material. 
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1. Introduction 

In order to advance our understanding of the ambient conditions 
of star formation, observations of large areas of known mas¬ 
sive Galactic star-forming regions have been done over a wide 
range of wavelengths. Observations of low- and mid-/ transi¬ 
tions of 12 CO towards several massive star-forming regions have 
shown that warm and dense gas is usually confined to narrow 
(< 1 pc) zones close to the ionization front (e.g., |Harris et ak 


A vast amount of data also shows that molecular clouds are 
clumpy over a wide range of scales, from giant molecular clouds 
containing smaller subclouds, to small dense cores within the 
subclouds. The dumpiness of molecular clouds is relevant for 
the understanding of fragmentation processes that can lead to 
compact condensations that may collapse to form stars (e.g., 
Carr||l987(|Banerjee et al.|2004l |Klessen et al.||2005[|Hocuk & 


1987 


2004 


Graf et al.||1993l [Y amamot o et afl|2001 [ |Kramer et al. 


|2008[ Perez-Beaupuits et al.|2010 ). Although slow shocks 


Spa ans|2010||Clark et al.|2011[[Federrath & Klessen|2012| and 


and cloud-cloud collisions can be an important source of heating 
in high velocity wing objects like Orion, W51, and W49 ( |Jaffe| 
narrow mid-/ 12 CO lines, as well as the parame- 


references therein). Several efforts have been made over the past 
years to identify clumps, their size, line width, and mass, among 
other parameters. The complex line profiles observed in optically 
thin lines (e.g., [C i], CS, C 18 0, 13 CO) and their velocity-channel 
maps, are indicative of the clumpy structure of molecular clouds 


et al. 


ters needed to explain the CO observations, favor photoelectric 

Carr 1987; Loren 1989; 

Stutzki & Glisten 1990 

Hobson 1992; 

heating of the warm gas located beyond the H n region driven 

Kramer et al. 19981 2004 

| Perez-Beaupuits et al. 

201 Of - 


by the UV radiation field emerging from an ionizing source, the 
so-called photon-dominated region (PDR). 


Send offprint requests to: J.R Perez-Beaupuits 


Inhomogeneous and clumpy clouds, as well as a partial 
face-on illumination, in star-forming regions like M17 SW, 
NGC3603, S140, the Orion Molecular Cloud, and the NGC 7023 
Nebula, are thought to produce extended emission of the atomic 
lines [C i] and [C n], and suppress the stratification in [C n], [C i], 
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and CO expected from standard 1-D steady-state PDR models, 


which is not observed in several sources (e.g., 

Keene et al. 1985, 

Genzel et al. 1988; Stutzki et al. 1988; Spaans & van Dishoeck 

1997 

Gerin & Phillips 1998; Yamamoto et al. 2001][ Schneider 

et al. 

[2002; Rollig et al. 2011; Mookerjea et al. 2003, Perez- 

Beau] 

puits et al. 2010 ). 


Massive star-forming regions like the Omega Nebula M17, 
with an edge-on view (particularly in its southwest region), are 
ideal sources to study the clumpy structure of molecular clouds, 
as well as the chemical and thermodynamic effects of the nearby 
ionizing sources. The southwest region of M17 (Ml7 SW) con¬ 
centrates molecular material in a clumpy structure. Models based 
on far-IR and submillimeter observations (Stutz kTet al.| [l988; 
Mei xner et al.||1992| ) suggest that the distribution and intensity 
of the emissions observed in the M17 SW complex, can be ex¬ 
plained with high density (ft(H 2 ) ~ 5 x 10 5 cm -3 ) clumps em¬ 
bedded in an interclump medium (n( H 2 ) ~ 3 x 10 3 cm -3 ) and 
surrounded by a diffuse halo (n( H 2 ) ~ 300 cm -3 ). 

The central cluster of more than 100 stars that illuminates 
M17 SW is NGC 6618 (e.g., |Lada et~ahl[T99T[ |Hanson et al 


1997|) . The two components of the massive binary CEN1 ( [Klein- 
1973[ |Chini et aL| [T980) are part of the central cluster 


mann 


NGC 6618 and are separated by ~ 178. This source, originally 
classified as a double O or early B system by Kleinmann (1973), 
is actually composed of two 04 visual binary stars, named CEN 
la (NE component) and CEN lb (SW component), and it ap¬ 
pears to be the dominant source of photo-ionization in the whole 
Ml7 region ( [Hoffmeister et al.|2008| ). 

Recent SOFIA/GREAT observations of the velocity - 
resolved [C n] spectra showed that a large fraction (>60%) of the 
[C n] emission, observed at the lower (<10 kms -1 ) and higher 
(>24 kms -1 ) velocity channels, is not associated with the star¬ 
forming material (denser and colder gas) traced by species like 
CO and [C i], which has an average line width of 5 to 10 k ms 


centered at Vlsr = 20 kms (Perez-Beaupuits et al 


2012 


2013). Only the central narrow (1 kms L ) channel maps of the 


velocity-resolved [C n] spectra show a spatial association with 
other gas tracers (e.g., [C i] and 12 CO). The broader velocity 
range covered by the [C n] line with respect to the [C i] and 12 CO 
has to be associated with additional material, either lower den¬ 
sity clumps or more diffuse, possibly ablated material, resulting 
in additional layers of ionized carbon gas within the telescope 
beam. The [C n] emission have been found to extend at least | 
in the sky ( [Russell et al.[|1981| , and ~15 pc into the Ml7 SW 


molecular region (Stutzki et al. 


1988| ). The spatial distribution 


of the [C n] emission (and abundance) in the southern region of 
Ml7 SW does not follow theoretical predictions of stratified or 
clumpy PDR models (P erez- Beaupuits et al J2012[ ). 

In earlier works, high resolution maps of high- and mid-/ 
CO lines, the 3 P 2 —> 3 P i fine-structure transition of [C i], and 
the [C n] 158 pm emission, have been reported (Perez-Beaupuits 
et al.|20lO||2012| ). In this study we present a new high resolution 
map of the ' P\ —> 3 Po fine-structure transition of [C i], as well 
as maps of the / = 1 —> 0 and J = 2 —> 1 transitions of 12 CO and 
its isotopologues. In contrast to [C n] 158 pm (and [O i] 63 pm, 
not included in the present data set), PDR models predict that 
the intensity of the [C i] fine structure lines do not have a strong 
dependence on UV energy density (e.g., [Hollenbach & Tielensj 
|1999| . Therefore, in a clumpy cloud irradiated by UV photons, 
the intensity of the [C i] emission is expected to be proportional 
to the number of p hotodissociation surfaces of clumps along the 
line o f sight (e.g., |Spaans|1996[ [Howe et al.|2000| [Kramer et al.| 
2004). Since several velocity components along the line of sight 


can be found in a clumpy medium, we present our analysis and 
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discussions of the new results based on velocity channel maps, 
showing the temperatures of the lines integrated over a narrow 
1 kms -1 channel width. From them we estimate the excitation 
temperature and column density of [C i], as well as the column 
density of [C n] and the gas mass not associated with the star¬ 
forming material traced by [C i] and the CO isotopologues. 

The organization of this article is as follows. In Sect. 2 we 
describe the observations. The maps of the observed lines are 
presented in Sect. 3. The excitation temperature and column den¬ 
sities, as well as mass estimates, are presented in Sect. 4. In Sect. 
5 we estimate the [C n] emission not associated with other gas 
tracers. The conclusions and final remarks are presented in Sect. 
6. 


2. Observations 


2.1. The APEX data 

We used the higher frequency band of the dual channel receiver 
FLASH ( [Heyminck et aT7[ pD06| hereafter FLASH-460) on the 
Atacama Pathfinder Experiment (APEXF Gusten et al.| [2006) 
during October 2009 to map the 3 Pi —> 609 pm (hereafter 

1 —> 0) fine-structure transition of [C i] at 492.161 GHz. The ob¬ 
served region covers about 672 x 772 (4.1 pc x 4.7 pc) compared 
to the 5'.3 x 477 (3.4 pc x 3.0 pc) area previously mapped for 
[C i ]JP 2 -» 3 Pi 370 pm (hereafter / = 2 -> 1) with CHAMP + 
(Perez-Beaupuits et al. 2010[ ). The [C i] J = 1 —> 0 was ob¬ 
served in on-the-fly (OTF) slews in RA (~ 360 arcsec long). 
Because the beam size of APEX at 492 GHz is about 1277, the 
subsequent scans in Declination and RA were spaced 6" apart. 

The total power mode was used for the observations, nodding 
the antenna prior to each OTF and raster slew to an off-source po¬ 
sition (180",0"), east of the star SAO 161357. This is used as the 
reference position (A a = 0, AS = 0) in the maps and through¬ 
out the paper, with RA(J2000)=18:20:27.64 and Dec(J2000)=- 
16:12:00.90. The OFF position at 180" was determined to be 
clean, even in the [C i] J = 2 —> 1 and mid-/ 12 CO lines. The 
reference for continuum pointing was Sgr B2(N) and the point¬ 
ing accuracy was better than 3" for all the maps. The data were 
processed during the observations with the APEX real-time cal¬ 
ibration software ( [Muders et al.|2006] ), assuming equal gains for 
the signal and image sidebands. 

A fast Fourier transform spectrometer (FFTS), providing 1.5 
GHz bandwidth and 2048 channels ( [Klein et al.|2012| ), was used 
for the [C i] / = 1 —> 0 map. The on-source integration time per 
dump was 1 second for the OTF map of [C i] / = 1 —> 0, and the 
average DSB system noise temperature of the FLASH-460 was 
about 810 K. 

Observations toward Jupiter were performed during October 
2009 to estimate the beam coupling efficiency (r] c « 0.59) of 
the FLASH-460, assuming a brightness temperature of 158 K 
for this planet at 492 GHz, as interpolated from data reported in 
[Griffin et al.[ ( 1986 ). 


2.2. The IRAM 30m data 


We used four frequency setups of the broadband EMIR receivers 
( Carter et al.||2012| ) at IRAM 30m to map a similar area to the 
one mapped in [C i] with the APEX telescope. The 32 GHz sig¬ 
nal bandwidth provided by the IF channels were used for the 


1 This publication is based on data acquired with the Atacama 
Pathfinder Experiment (APEX). APEX is a collaboration between the 
Max-Planck-Institut fiir Radioastronomie, the European Southern Ob¬ 
servatory, and the Onsala Space Observatory 
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In order to ensure atmospheric stability, we used the 
same nearby off-source reference position (345",-230") as for 
the [C n] map (Perez-Beaupuits et al] |2012| ) obtained with 
the Germa n REceiver for A stronomy at Terahertz frequencies 
(GREAI^J |Heyminck et al. 2012 ) on board the Stratospheric 
Observatory For Infrared Astronomy (SOFIA). From previous 
APEX observations (not reported here) of the 12 CO J - 3 —> 2 
and / = 4 —> 3, we know that the reference position (345",- 
230") is not free of CO emission. Hence, for all the EMIR 
frequency setups, we first did a deep observation at the po¬ 
sition (345",-230") against the reference that is even farther 
away at (3600", -1800") which, in turn, was verified to be CO 
emission-free against the offset position at (4600",-2800"). Then 
we added the flux from the reference position (with about two or¬ 
ders of magnitude lower rms than the OTF spectra) back into the 
spectra from the OTF maps. 

Beam coupling efficiencies for each individual line were ob¬ 
tained from interpolation of the values given in the online table 
of the IRAM 30m efficiencie^] With these beam coupling effi¬ 
ciencies, and a forward efficiency pf (0.95 for EMIR090, 0.94 
for EMIR230, and 0.93 for EMIR150), we converted all data to 
the main beam brightness temperature scale, T B = 

The reduction of these calibrated data, as well as the maps shown 
throughout the paper, were done using the GILDAS^ package 
CFASS90. 


2.3. SOFIA [C n] observations 


The [C n] 158 pm map was already reported in Perez-Beaupuits 
et al.| ( |2012| ), where a detailed explanation of the calibration and 
the OTF maps was given. Here we use the same data as in the 
previous paper, convolved with a larger beam, and re-sampled 
the spectra to a 1 km s -1 channel width, as explained in Sect. [5] 
The SOFIA data is publicly available in the section “Data 
Archive & Retrieval” of the SOFIA Data Cycle Systerrj^] All the 
data presented in this work will be available as FITS files in the 
Strasbourg astronomical Data Center (CDS^J. 


Fig. 1 . Top - Color map of the velocity integrated (in the range 0- 
40 km s' 1 ) intensity of [C i] / = 1 -» 0 in M17 SW. The peak emission 
is 240 Kkms -1 . The contour levels are 25%, 50%, 75%, and 90% of the 
peak emission. Bottom - Color map of the integrated intensity of [C i] 
/ = 2 —» 1 (from Perez-Beaupuits et al. 2010) convolved to the beam 
size (~ 1277) of the [C i] / = 1 —» 0 line, with a peak emission of 
280 Kkm s -1 . The contour levels are as described above. The reference 
position (A a = 0, AS = 0) is as in Fig. [2] 


receiver bands E090 (3mm) and E230 (1.3mm), covering each 
sideband with 8 GHz bandwidth in single polarization. These 
setups allowed us to fully map all the CO isotopologues (in addi¬ 
tion to many other molecules) in their / = 1 —» 0 and J - 2 —> 1 
transitions. The beamwidths (FWHM) for the / = 1 —> 0 tran¬ 
sitions of 12 CO, 13 CO, C 18 0, and C 17 0 are 22".6, 23".7, 23".8, 
and 23".2, respectively. We also detected, by serendipity, the hy¬ 
drogen recombination lines H39 a, H40 a, and H41 a (2873) in 
the 3mm band. 

The total region mapped of about 360"x300" was covered 
with two long OTF maps of 360"x 160" (with an overlap of 20" 
between them) and slews in RA (~ 360 arcsec long) with steps 
of 4" in Dec. The off-source reference position was observed 
for 10 seconds every two OTF subscans (rows). The on-source 
integration time per dump was 0.5 second. 


3. Results 


3.1. The [C i] integrated intensity maps 


Figure [I] shows the maps of the intensity, integrated between 
0 km s^and 40 km s -1 , of [C i] / = 1 —> 0 (top) and J = 2 —> 1 


(bottom). Because the [C i] / = 2 —> 1 map was convolved to 
the larger beam size (1277) of the / = 1 —> 0 transition, its peak 
integrated intensity is ~ 20 K km s -1 l ower than the peak value 
previously reported in |Perez-Beaupuits et al.[ ( [20T0| . The peak 
integrated intensities of the maps shown here are 240 Kkms -1 
and 280 Kkms -1 for the [C i] / = 1 —> 0 and J - 2 —> 1 lines, 
respectively. 

These lines follow a similar spatial distribution and their re¬ 
spective peaks are located at about the offset position A a = 
-120", A S = 30", approximately 0.88 pc (~ 80" at RA. 90°) 
from the ridge. They both present extended emission, unlike a 
theoretically expected stratified PDR . A spatial associati on be¬ 
tween 13 CO and [C i] was found by 


Keene et al. 


(1985), 


on a 


2 GREAT is a development by the MPI ftir Radioastronomie and the 
KOSMA/ Universitat zu Koln, in cooperation with the MPI ftir Sonnen- 
systemforschung and the DLR Institut ftir Planetenforschung. 

3 http://www.iram.es/IRAMES/mainWiki/Iram30mEfRciencies 

4 http://www.iram.fr/IRAMFR/GILDAS 

5 https://dcs.sofia.usra.edu/ 

6 http://cdsweb.u-strasbg.fr/ 
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Fig. 2. Color maps of the velocity integrated intensity (in the range 0-40 km s l ) of the J = 1 —» 0 and / = 2 —» 1 transitions of 12 CO (Top panels) 
and 13 CO (Bottom panels) in M17 SW. The contour levels are 25%, 50% (thick contour), 75%, and 90% of the peak emission. All maps have 
been convolved to the larger beam (24") of the C 18 0 / = 1 —» 0 line. The reference position (A a = 0, A 8 - 0), marked with a cross, corresponds 
to the SAO star 161357 at RA(J2000)=18:20:27.65 and Dec(J2000)=-16:12:00.91. The ultracompact H n region M17-UC1 and four H 2 0 masers 
( [Johnson et al.|1998) are marked by the black circle and plus symbols, respectively. 


scale size of 3' (~ 2 pc). We confirm these results with higher 
spatial resolution. We also add that, when looking at the overall 
distribution, the C 17 0 and C 18 0 emissions are more similar to 
the 13 CO than to the 12 CO emission and, hence, they also show 
a spatial association with the [C i] integrated emission. This spa¬ 
tial association is discussed further in the next section. 


3.2. The CO integrated intensity maps 

In order to compare them with the J = 1 —> 0 transitions, and to 
increase the S/N of the J = 2 —> 1 lines, all the maps (including 
the [C i] lines from Fig. [T]) were convolved to the larger beam 
size (24") of the C 18 O J = 1 —> 0 line for the analysis presented 
in the next sections. Maps of the velocity integrated intensity of 
the 12 CO emission, and its isotopologues 13 CO, C 17 0, and C 18 0 
lines, are shown in Figs.[2]an d0 

The 12 CO integrated emission is more extended, and its bulk 
emission does not resemble that of its isotope lines. In order 
to verify this, we use the scatter plots (Fig. [4} of the velocity- 
integrated intensity of these tracers, and the corresponding cor¬ 
relation coefficient described in Appendix [C] The scatter plots 
deviate from the theoretical straight line expected for well- 
correlated maps. The 13 CO/ 12 CO plots show a well-separated 
optically thin (low intensities) and an optically thick (high in¬ 
tensities) branch. It is interesting to note that the optically thick 
branch still has a relatively good correlation. This shows the in¬ 
tegrated intensity growth via line broadening at lower densities, 
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fully in line with Larson’s law; i.e., the cloud size is inversely 
proportional to its density, and the velocity dispersion (or line 
broadening) is directly proportional to the cloud size, hence, the 
lower the density, the larger the cloud and the broader the lines. 

We note, however, that even though the intensities of the 
12 CO and 13 CO maps are very scattered (at the higher values) and 
show two branches with different slopes, the correlation coeffi¬ 
cient is still relatively high, r xy = 0.94 - 0.96. These correlation 
coefficients are similar to those found between the / = 1 —> 0 
and J = 2 —> 1 lines of 12 CO and [C i], which are much better 
correlated, as shown in Fig. |C.1| A similar case is found when 
comparing the [C i] 370 pm line with the / = 2 —> 1 transi¬ 
tions of the 13 CO and C 18 0 lines (bottom panel in Fig. [4]). These 
scatter plots show values that are less scattered (supporting the 
similarity between the [C 1 ] and the CO isotopologue lines men¬ 
tioned above), while the correlation coefficient is similar to (or 
lower than) those found for maps with less similar spatial distri¬ 
bution (e.g., 12 CO and 13 CO J = 1 —> 0). This means the cor¬ 
relation coefficient is not robust enough to discriminate between 
“well” correlated and “not so well” correlated maps and, hence, 
it must be used with caution. Another way of using this statistical 
measure is to define a threshold to consider only regions of the 
map with high intensities, to compute the correlation coefficient 
as described in Sect. 15.11 

Figure [5] shows the overlay between the two J = 1 —> 
0 and 2 —> 1 lowest energy transitions of 12 CO and the 
[C 11 ] 158 pm line, both in the total velocity-integrated inten¬ 
sity and the 1 km s _1 channel maps. We note that the velocity- 
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Fig. 3. Color maps of the velocity integrated intensity (in the range 0-40 km s l ) of the J = 1 —» 0 and J - 2 —> 1 transitions of C 17 0 (Top panels ) 
and C 18 0 (Bottom panels ) in M17 SW. Only the C 17 0 7 = 1 —> 0 line was integrated between 14 and 28 kms -1 because it is fainter and noisier 
than the other lines. The contour levels are 25%, 50% (thick contour), 75%, and 90% of the peak emission. All maps have been convolved to the 
larger beam (24") of the C 18 0 / = 1 —» 0 line. The reference position (A a - 0, AS = 0) and symbols are as in Fig. [2] 


integrated intensity of [C ii] does not follow the spatial distri¬ 
bution of the 12 CO emission. However, some spatial association 
between these lines can be seen in the channel maps, but only at 
the central 16-24 km s -1 velocity channels, where the bulk of the 
12 CO emission is found. 


New results in Ml7 SW have shown that a large fraction 
of the [C n] emission is not associated with other species (e.g., 
12 CO and [C i]) tracing the star-forming material when analyzed 
at narrow velocity channels ( [Perez-Beaupuits et al.|2012[[2013| ). 
Hence, line integrated intensity maps have to be interpreted with 
great care, as a smaller or larger part of the line integrated emis¬ 
sion, as in the case of [C n] here, may result in strong veloc¬ 
ity components in one line, which are barely traceable in an¬ 
other line, like C 18 0 or 13 CO here, due to different physical ori¬ 
gin and/or very different excitation conditions. This obviously 
also strongly affects the interpretation of line-integrated inten¬ 
sity ratios between different tracers. In fact, we do not find 
strong spatial associations when comparing the velocity inte¬ 
grated [C n] 158 pm map from Perez-Beaupuits et al. (2012, 
their Fig. 1) with the new CO and [C i] maps from Figs.[2]an d0 
Therefore, in the following sections we present our new analysis 
and discussions based on velocity channel maps, showing the in¬ 
tensities of the lines integrated over a narrow (1 km s -1 ) channel 
width. 


4. Excitation and column density of [C i] 

The critical density (n cr ~ 10 3 cm -3 for collisions with o-/p- 
H 2 at 100 K; from the LAMDAp] database, |Schoier et al.|2 005 ) 
and upper-level energy (E u ^ 24 K for J = 1 —> 0, and E u « 
62 K for J = 2 —> 1) of [C i] enable us to trace the diffuse 
(n( H 2 ) < 10 3 cm -3 ) ISM and estimate its temperature. We note, 
however, that this does not imply that [C i] traces only diffuse 
gas, but also denser gas given the overall spatial association of 
the velocity-integrated intensities observed between [C i] and the 
CO isotopologues, as described in section |3H] 

We first estimate the excitation temperature of [C i] from the 
ratio between the two transitions, assuming optically thin emis¬ 
sion. From this excitation temperature, the column density of 
[C i] can be estimated as well, assuming optically thin emission 
and LTE conditions. Then we compare the results with a non- 
LTE estimate at representative offset positions. This is done in 
section l4~2l 

4.1. Excitation temperature of [C i] 

We computed the ratio R = 7([C i] 369/im)//([C i] 609yum) 
between the intensities of the [C i] lines integrated in narrow 
(1 km s -1 ) velocity channels, for each channel map in the veloc¬ 
ity range between 0 km s -1 and 40 km s -1 . We find observed val¬ 
ues of 1<R<2, at the central velocity channels (16-24 kms -1 ). 
Ratios lower than unity (which means subthermal excitation) 

7 http://www. strw. leidenuniv.nl/~moldata/ 
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As mentioned in Sect. 3.1 the total integrated intensity of 
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Fig. 4. Scatter plots and correlation coefficients r xy (from Eq. |C.1| > be¬ 
tween the velocity-integrated intensity of 12 CO and 13 CO J = 1 —» 0 
(top-left), and the 12 CO / = 2 —> 1 compared to the / = 2 —> 1 lines 
of 13 CO ( top-right ), C 17 0 (imiddle-left ), and C 18 0 (middle-right). The 
bottom panel shows the scatter plot and correlation coefficient between 
the [C I] 370 pm and the / = 2 —> 1 transitions of 13 CO (bottom-left) 
and C 18 0 (bottom-right). 


are also observed in our maps, but mostly at velocity channels 
<16 kms -1 and >24 kms -1 , where the intensity of both lines is 
below 30% of their respective peak channel intensities, which is 
about the noise level of the fainter [C i] 609yum line. 

The observed ratios match the values expected in a PDR en¬ 
vironment with low density 10 3 cm -3 ) and relatively low radi¬ 
ation fields (Go £ 10 2 , with Go in units of the ambient interstel¬ 
lar radiation field, 1.2 x 10 -4 ergs s -1 cm -1 sr -1 , Habing 1968 ), 
as shown by Meijerin k et al.|p507] their Fig. 3). These densi¬ 
ties are also in agreement with previous estimates by |Mei xner 


et al. 


(2007 their Fig. 3) at higher densities (~ 10 4 cm 3 ) and 


even lower radiation fields (Go ^10), which are naturally atten¬ 
uated by the larger column densities of denser clumps. 
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the [C i] emission resembles that of the optically thin isotope 
CO lines, rather than that of the 12 CO lines. The isotope lines 
have critical densities of n crit ~ 10 4 cm -3 (at T K 100-200 K) and 
trace the compact emission from the denser molecular clouds. 
Because of the self-pumping due to large optical depths result¬ 
ing from the high abundance, the 12 CO lines trace not only 
the denser molecular clouds, but also the more diffuse and ex¬ 
tended emission beyond (i.e., east from) the ionization front. 
Hence, contrary to the general picture proposed by |Meixner et al.| 
(19921, we favor a scenario where part of the [C i] emission 


|et a 1. |(T992] ), who concluded that the [C n], [C i], and low -J CO 
lines emerge from an inter-clump medium, while the extended 
[C n] and [C i] emission emerges from a halo gas surrounding 
the clump and inter-clump material. However, the observed ra¬ 
tios are not exclusive of these ambient conditions, since they can 
also be foun d by extrapolating the values given by Meij erink] 


emerges from an inter-clump medium rather than from a more 
diffuse halo gas, and the other part is associated to the denser 
clumps traced by the isotope CO lines. 

Following the formalism by Schneid er et al.| ([2003 , their Ap¬ 
pendix A), which is valid in the optically thin limit and when 
both lines have similar optical depths, we can estimate the ex¬ 
citation temperature T ex of [C i] from the ratio R as T ex = 
38.8/7/2(2.11//?) K. Figure [6] shows the estimated T ex of [C i] at 
the same region mapped in [C n] for the velocity channels from 
9 kms -1 to 29 kms -1 where the emission of both [C i] lines is 
larger than the 3 cr level detection. The [C i] maps were first con¬ 
volved to the larger beam (24") of the C 18 O J - 1 —> 0 map in 
order to increase the S/N. 

The excitation temperature ranges between ~40 K and 
~100 K in the inner molecular region (i.e., southwest from the 
ionization front). This result is in agreement with a previous LTE 
estimate of the excitation temperature from mid -J 12 CO (Perez- 
Beaupuits et al. 1)2010 ), earlier estimates from [C i] / = 2 —> 1 
observations (e.g., |Genzel et al.||1988| ) that indicated the [C i] 
emission arises from gas with kinetic temperature of about 50 K, 
and from a multi-line NH 3 study ( [Glisten & Fiebig|1988| ) which 
showed different coexisting gas phases with kinetic temperatures 
between 30 K and few 100 K, and up to about 275 K in the re¬ 
gion traced by the VLA continuum arc (the northern ionization 
front). Even higher excitation temperatures are found at sparse 
locations northeast of the ionization front, and along the east¬ 
ern edge of M17 SW at the channel bins (20-24 kms -1 ) close 
to the ionization front. Pointing offsets between the CHAMP + 
and FLASH observations (done at different observing periods), 
as well as differential couplings of the two respective beams just 
at the edge of the ridge, can mimic such temperature gradients. 
These effects could account for up to 30% in the line ratios, 
and they cannot be discarded. Hence, we limit the color scale 
of Fig.[S]to 200 K, which is a more reliable upper limit of the ex¬ 
citation temperature in the region mapped. Ratios >2, found at 
few positions between the ionization front and the ionizing stars, 
lead to negative excitation temperatures in the optically thin and 
LTE approximation. 


4.2. Optical depths and column density of [C i] 

From the excitation temperature and the peak intensity of the 
[C i] / = 1 —> 0 and J = 2 —> 1 lines, the optically thin ap¬ 
proximation also allows us to estimate the optical depths of both 
lines. Knowing the excitation temperature and the optical depth 
of the / = 1 —> 0 line, the column density A([C_i]) can be com¬ 
puted as well. For detailed formulae see [Frerking et al.| ( |1989| ) 
and Schneider et al. (2003, their Appendix A). We also assume 
that the sources fully cover our beam in the emitting regions (i.e., 
we use a beam filling factor of unity). Therefore, the quantities 
reported here correspond to beam (24") averaged values. 

Figures |A.1| and |A.2| show the channel maps of the optical 
depths estimated for the [C i] / = 1 —> 0 and / = 2 —> 1 
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Fig. 5. Top - Velocity-integrated intensity maps of 12 CO / = 1 —» 0 (gray), [C n] 158 pm (red contour), and 12 CO / = 2 -> 1 (green contour). The 
contour lines_(from thin t o thick) are 50%, 75%, and 90% of the respective peak emissions. The stars indicate the O and B ionizing stars ( jBeetz| 
|et al.|1976||Hanson et al.|1997] >. The reference position (Acr = 0, AS = 0), marked with a cross, is as in Fig. [2] Bottom - Velocity channel maps (at 
1 kms _i width) of the same lines as above. Contours are 20%, 40%, 60%, 80%, and 100% of the respective peak emissions. All maps have been 
convolved with the largest beam of 24" corresponding to the C 18 0 / = 1 —» 0 map. 


lines, respectively. Optical depths r < 1 are observed in most 
of the regions mapped at all the velocity channels. Except in 
a small region around the offset position (-130", -10") of the 
T 2 —>i line, at the central velocity channels (19-20 km s -1 and 20- 
21 km s _1 ), where the bulk of the [C i] emission is found. 


A non-LTE exc itation analysis using the Radex code ( [van 
der Tak et al.|2007| ), was used to test the optically thin assump 

ti rvvi 1 n 4 11 r/\ /-4-n + n+iT 7/-^ rvvin i-m tVi n 1 O O M 1 / n 1 n-M 


tion in two representative positions in the 19-20 kms chan¬ 
nel map. First at offset position (-130", -10") where T 2 ->i > 1 


(cf. Fig. B.l| ), and then at the offset position (-130", -70") (cf. 
Fig. |B.2 ) where T 2 ->i < 1, as seen in Fig. [7] The [C i] line 
ratios and intensities observed at these positions can be repro¬ 
duced with densities larger than 10 3 cm -3 (the critical density of 
[C i] / = 1 —> 0 at T k between 100 K and 200 K) for kinetic 
temperatures below 500 K, and column densities per line width 
N/AV = (4 - 7) x 10 17 cm -2 km -1 s, similar to the column densi¬ 
ties obtained with the FTE method. Considering a line width of 
~ 10 km s _1 , the column densities we obtained are consistent with 
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Fig. 6. Velocity channel maps (at 1 kms -1 width) of the [C i] excitation temperature (color map in K) in M17 SW, estimated from the R = 
7([C I] 369ji/m)//([C I] 609yum) line intensity ratio, and assuming LTE conditions and optically thin emission. Only velocity channels and pixels 
with emission larger than 3cr in both [C I] lines are shown. 



100 0 -100 -3GK* 

A* [ - ] 



A* [ " 1 


Fig. 7. Optical depths of the [C i] / = 1 —» 0 (top) and / = 2 —» 1 ( bot¬ 
tom )i lines i n the 19-20 kms -1 velocity channel map toward M17 SW 
(from Figs. |A.l| and |A.2| ), estimated from the excitation temperature of 
[C i], and assuming LTE conditions. 


the peak column density and a moderate optical depth of about 
2.5+0.7 derived by |Genzel et al. ( 1988] ) from the total velocity- 
integrated intensity of [C i] J = 2 —> 1, and assuming a non¬ 
linear relation between the intensity of [C i] and the intensity of 
C 18 0. 


Our observed [C i] intensities can also be reproduced with 
densities above ~3xl0 3 cm -3 (the critical density of [C i] J = 
2 —> 1 at Tk =100-200 K) and temperatures below 300 K. The 
[C i] / = 1 —> 0 line is close to thermalized (T ex « 7\) at both 
positions, for T k < 500 K and densities n( H 2 ) = 10 3 - 10 4 cm -3 , 
while r is not much smaller than unity in both positions, just 
as in the LTE results. Although the optical depths of both [C i] 
lines are just marginally thin, they are very similar. Hence, we 
can still apply the optically thin approximation in all the regions 
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of interest, at least for which is the opacity used to estimate 
the column density of [C i] following [Schneider et al.| ( 2003 


their Eq. A.8). Since the optical depth is not much smaller than 
unity, we used the correction factor r([C i])/(l - exp(-r ([C i]))) 
to compute the column density of [C i]. 

The velocity channel maps of the column density N([C i]) 

The values of N([C i]) range 


(cm 2 ) are shown in Fig. 
between 10 15 cm -2 and 


A.3 

W 7 ” 


cm -2 throughout the whole re¬ 
gion mapped and among all the velocity channels. However, the 
bulk of the [C i] emission corresponds to column densities above 
10 16 cm -2 in all the velocity channels. Column densities up to 
10 17 cm -2 are reached only in the central velocity channels, be¬ 
tween 17 km s -1 and 22 km s -1 , which correspond to the regions 
where the integrated intensity maps (Fig. [T} show a strong [C i] 
emission (> 50% of the peak). 


5. Emission of [C n] not associated with other gas 
tracers 

In Perez-Beaupuits et al. (2013) we showed overlays of our pre¬ 
vious 12 CO J - 2 


1 map and the optical depth of HI from 
Brogan & Troland ( 2001), over the [C n] emission, also in chan¬ 
nel maps of 1 km s 1 widt 


L maps or i km s x width. We showed that only at intermediate 
(10-24 kms -1 ) velocities the [C n] emission presented strong 
spatial association with other molecular gas tracers. A strong 
spatial association is identified when the spatial distribution of 
a particular [C n] channel map is very similar, or adjacent, to 
that of another gas tracer (e.g., 12 CO J - 2 —> 1 or [C i]). On 
the other hand, at lower (<10 kms -1 ) and higher (>24 kms -1 ) 
velocity channels, the [C n] emission is mostly not associated 
with the other tracers of diffuse and dense gas. We note that “not 
associated” in this sense does not mean that we deal with physi¬ 
cally completely independent material. The overlay between the 
12 CO 2 —> 1 and [C n] velocity channel maps in Fig. [5] shows 
that in the outer velocity range the [C n] emission often shows 
halos and diffuse extensions around the denser clumps and fila¬ 
ments identifiable in the 12 CO channel maps, suggesting that the 
[C n] emission traces gas that has been ablating off the clump or 
filament surfaces; this [C n] emitting gas, however, is not visible 
in 12 CO, despite this association. 
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Fig. 8. Velocity channel map (integrated in 1 kms -1 ) of the [C n] 158yi/m emission (gray background), with overlays of the 3 P i —> 3 P 0 hne- 
structure transition of [C i] at 609yum (contours). The ionized and neutral carbon emissions are well correlated only at the intermediate velocities 
(16-24 kms -1 ). 


This can also be seen in the line shapes of the [C n], 
[C i] 609 /d m, C 18 0 J = 2 —> 1, and 12 CO 7=1-0 spectra 
shown in Fig.[lT]for different offset positions along a strip line at 
position angle (P.A.)=63°. The pointings of the spectra from the 
different tracer coincide within 2", i.e., sufficiently close con¬ 
sidering the smeared 24" beam resolution. We note that most 
of the lower and higher velocity channels of the [C n] line are 
not associated with any of the other lines, while the C 18 0 and 
[C i] lines are highly associated. The velocity channel maps of 
[C n] with overlays of the new maps of the [CO 609 /dm and 
12 CO / = 1 —> 0 lines, are shown in Figs. [8] and [9] respectively. 


The spatial distribution of each velocity channel in the 
[C i] 609 /tm emission follows a similar pattern to the maps pre¬ 
sented in |Perez-Beaupuits et al.| ( [2013| , although with a narrower 
(14-24 km s _1 ) velocity range where strong association with the 
[C n] emission is observed. 


If the [C n] emission was associated in all velocity channels 
with the more diffuse molecular gas (between 30 and 300 cm -3 , 
independent of temperature) traced by the 12 CO / = 1 —> 0 
line, the spatial association between [C n] and 12 CO / = 1 —> 0 
would be expected to be stronger than that between [C n] and 
12 CO J = 2 —> 1. However, like the 12 CO J = 2 —> 1 and 
[C i] lines, the 12 CO / = 1 —> 0 line shows a strong associa¬ 
tion with [C n] only in the central 10-24 kms -1 components (cf. 
Fig. [9|. This is another confirmation that the lower (<10 km s -1 ) 


and higher (>24 km s ] ) velocity channels of the [C 11 ] emission 
are not strongly associated with the bulk of the molecular gas. 

5.1. Spatial correlation between the star-forming material 
and the [C 11 ] emission 

We quantify the spatial association observed at each velocity 
channel between the [C 11 ] emission and other gas tracers (e.g., 
12 CO J - 1 - 0, [C 1 ] 609 /rm) according to the procedure de¬ 
scribed in Appendix[D| We use 10% of the global (i.e., among all 
the channels) peak emission of each tracer as a threshold to iden¬ 
tify the pixels with significant emission to be used. This method 
allows us to estimate where in the region and by what fraction of 
the region the [C 11 ] emission is associated, in the spatial distri¬ 
bution of each velocity channel map, with other gas tracers (cf. 
first paragraph in Sect. [5]). 

The correlation coefficient described in Appendix[C] and his¬ 
tograms showing the velocity distribution of the spatial associa¬ 
tion between the [C 11 ] emission and that of [C 1 ] 60 9 yum, 12 CO, 
13 CO, and C 18 0 / = 2 —> 1 are shown in Fig.|lO[ In this case 
the correlation coefficient r xy was computed using only the pix¬ 
els with emission larger than 1/3 of the global peak emission, in 
order to consider only the optically thick branch of the CO scat¬ 
ter plots shown in Fig. [4] and to exclude the [C 1 ] scatter at low 
intensities. Figure [TO] then shows that the [C 11 ] emission is asso¬ 
ciated with, and highly correlated with (r xy >0.6) other gas trac- 
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Fig. 9. Velocity channel map (integrated in 1 kms -1 ) of the [C n] 158//m emission (gray background), with overlays of the 12 CO / = 1 —> 0 
transition (contours). The molecular gas shows better association with the ionized carbon at more extended intermediate velocities (14-27 km s -1 ) 
than the neutral carbon (Fig. [8]). 


ers, mostly at the central velocity channels between 15 kms -1 
and 23 km s -1 . The fraction of the mapped region where the [C n] 
emission is associated with the [C i] 609 pm line is 30%-55% 
in the velocity range mentioned above, while it reaches 40%- 
80% with 12 CO, 35%-50% with 13 CO, and only 20%-45% with 
C 18 0 J = 2 —> 1. 


The large range of velocity channels of [C n] emission not 
associated with other gas tracers is strong evidence of the in¬ 
ability of the total velocity-integrated [C n] line intensity (i.e., 
its total flux) to estimate its beam averaged abundance (i.e., col¬ 
umn density ratio in comparison to other tracers), and the cool¬ 
ing it provides to the molecular and star-forming gas traced by 
species like [C i] and CO. The assumption of the [C n] emission 
arising from the same spatial region as other tracers of diffuse 
and dense gas in the whole velocity range covered by the [C n] 
spectra, was the basis for many previous studies, according to the 
technology available at that time. In particular, the spectrometers 
on board of NASA’s Kuiper Airborne Observatory (KAO) had a 
spect ral resolution of 80- 175 kms -1 for the [C n] 158 pm line 
(e.g., 


Meixner et al. 


1992), while the [C n] spans only ~40 km s -1 
in our velocity-resolved spectra. Now we know that only ~20% 
of that resolved velocity range is associated with the star-forming 
material traced by [C i] and CO. Hence, using the total velocity- 
integrated [C n] intensity can be misleading, since it may include 
(as in the case of M17 SW) emission from gas that is not really 
associated with the [C i] or the CO emitting gas. Therefore, the 
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actual abundance of C + , and the cooling of the molecular gas 
due to [C n] emission associated with molecular gas in several 
Galactic and extra-galactic environments, may be overestimated. 

At the central velocity channels (e.g., 21-22 kms -1 ) the 
[C n] line shows dips at positions where optically thin lines (e.g., 
C 18 0 / = 2 —> 1) show their peaks, like the spectra in Fig, [ill 
at offset position (-60", -30"). This spectral feature is present at 
several positions around the previous one, indicating the pres¬ 
ence of a colder foreground layer, with significant optical depth 
to absorb the emission from a warmer background component 
present at an adjacent velocity. Taking the main beam temper¬ 
ature at the channel 22.5 kms -1 (where the C 18 0 / = 2 —> 1 
lines drops sharply) as the continuum level (7c ~ 65.3 K), and 
the temperature at the channel 21.5 kms -1 as the maximum ab¬ 
sorption depth of the [C n] line (T L ~ 47.4 K), we can estimate 
the optical depth of the absorbing layer as r = -ln( 1 - T^/Tq). 
Assuming that the absorbing layer completely covers the back¬ 
ground component and that all the foreground [C n] atoms are in 
the ground state, we can estimate the absorbing column density 
following Eq. (3) in |Neufeld et al.|p0l^ ) as 

f rdv = A » <g » /l3 jv(c + ) = 7.2x 10~ IS N(C + ) (cm- 2 ) km s -1 , (1) 
J 8 ngi 

where A u \ - 2.3 x 10 -6 s -1 is the spontaneous radiative decay 
rate, g u - 4 and gi = 2 are the degeneracies of the upper and 
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Fig. 10. Histograms showing the fraction of the [C n] 158 /dm emitting region correlated to each 1 kms 1 channel with other gas tracers. The 
corresponding correlation coefficient r xy is overlaid. 


lower states, and A = 157.741 jum (used in cm in Eq. [TJ is the 
transition wavelength of the [C ii] ground state. From this we 
estimate a foreground absorbing column density of N(C + ) « 2 x 
10 17 cm -2 . This is about a factor of two smaller than the column 
densities that we estimate for the background emitting [C n] gas 
(see Sect. 5.2). Since the absorbing layer affects only one or two 
of the central 1 kms -1 wide velocity channels, given the sharp 
profile of the C 18 0 line, it does not affect any of our conclusions. 

In the following section we estimate in more detail the [C n] 
column density and atomic hydrogen mass not associated with 
the star-forming material traced by the [C i] and C 18 0 lines by 
using an LTE approximation. 

5.2. Column density of [C n] and mass of dissociated 
material 

Because the spectra of [C n] and the other gas tracers have 
slightly different velocity resolutions, we have re-sampled all the 
spectra to 1 kms -1 resolution. 

We use the [C i] 609 /im line as a tracer of the diffuse 
(«(H 2 )crit ~ 400 cm -3 at T K -100-200 K) molecular gas, and 
the optically thin C 18 0/ = 2—>1 line as a tracer of the denser 
(ft(H 2 ) C rit ~ 9 x 10 3 cm -3 at T K -100-200 K) molecular gas in 
the star-forming material. 

For each [C i] and C 18 0 spectra of the map we find the chan¬ 
nel with the maximum intensity, and divide the corresponding 
[C n] channel by that maximum value to obtain the factors by 
which to multiply separately the [C i] and C 18 0 spectra. This 
produced scaled up [C i] and C 18 0 spectra that match the orig¬ 
inal [C n] spectra at their respective maximum channel intensi¬ 
ties. Then we subtract the scaled up [C i] line from the original 
[C n] spectra. The same is done independently with the C 18 0 


line to produce two residual [C n] spectra. Since we do not ob¬ 
serve any absorption line profile in our spectra, we consider all 
the channels with negative values as noise, and we set them to 
zero in order to avoid unwanted boosting of the [C n] emission. 
If the maximum intensity of [C i] or C 18 0 is higher than the in¬ 
tensity of the corresponding channel in the [C n] spectrum, then 
no scaling up is done. 

An example of this procedure, and the results for the spec¬ 
tra at the approximated offset position (-30",-15"), is shown in 
Fig. IT2] The shaded histogram corresponds to the residual spec¬ 
tra after the subtraction of the [C i] and C 18 0 lines. For each 
channel we took the minimum intensity between the two resid¬ 
ual [C n] spectra to produce a synthetic spectrum that represents 
the [C n] emission not associated with the star-forming material 
traced by [C i] and C 18 0. All these assumptions, in particular 
the scaling of the intensities, will lead to conservative numbers 
(lower limits). 

We can estimate the column density N(C + ) of the non- 
associated [C n] gas from the synthetic residual [C n] spectra, 
following the high-temperature FTE limit, which is valid for 
temperatures well above 91 K and high densities, 


N (C + ) : 


ri l I[ cii]6.3 x 10 


l20 cm 2 


( 2 ) 


obtained from a two-level system model as described in |Schnei^| 
|der et al.| ( |2003| their Eq. (A.5),) where r] c is the beam filling 
factor assumed to be unity since the [C n] emission in M17 SW 
is very extended, and /[cii] is the [C n] emission in units of erg 
cm -2 s -1 sr _1 . We note that with the two-level system expression, 
the estimated column density of [C n] increases if lower densities 
and/or lower temperatures are used. If we assume a gas tempera¬ 
ture of 250 K and a density of 10 4 cm -3 , as estimated for the re- 
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Fig. 11. Spectra of several lines observed at approximated (±2") offset 
positions along the strip line at P.A. 63° (A 6 = Aor/2). All the spectra 
have been resampled to a 1 km s -1 resolution and convolved with the 
largest beam of 24", corresponding to the C 18 0 / = 1 —> 0 map. 


gion east of the ionization front in M17 SW (i.e., the H i and H ii 
regions) from mid -J 12 CO line observations (Perez-Beaupuits et 
al. 2010), the exact expression for the two-level system presented 
by |Schneider et al.| ( |2003| their Eq. (A.4),) gives a 25% larger 
column density than the LTE approximation of Eq. Since 
the procedure described above subtracts the [C n] emission as¬ 
sociated with most of the dense and diffuse molecular gas, the 
residual emission should be dominated by collisional excitation 
from atomic hydrogen and free electrons. 


In order to analyze how sensitive the [C ii] column density 
is to an assumed temperature and density of the gas, we es¬ 
timate N(C + ) for a range of temperatures and densities using 
Eq. (A.4) from Schneider et al. (2003), considering a filling fac- 



Fig. 12. Residual [C n] 158 pm spectrum (dashed line and gray filled 
histogram) at offset position (-30",-15") after subtracting the scaled up 
(marked with (*)) [C i] 609 pm (top) and C 18 0 / = 2 —> 1 (bottom) 
lines from the original [C Ii] spectrum. All negative (noise) channels in 
the scaled up and residual spectra are set to zero. The residual [C n] 
spectrum is shifted by -1 K for clarity. 


tor pc of unity. We consider the critical densities n cr for free elec¬ 
trons (e~), atomic hydrogen (H i), and molecular hydrogen (Eh), 
which were computed for each temperature, according to the 
corresponding collisional deexcitation rate coefficients reported 
by |Barinovs et al.| ( |2005| ). Since the actual column density of 
[C ii] also depends on the [C n] intensity (which is arbitrary for 
this analysis), in Fig. [13] we show only the ratio with respect to 
the N(C + ) obtained using the temperature (250 K) and density 
(ft(El) = 10 4 cm -3 ) assumed above, to demonstrate the relative 
effect of using different ambient conditions. Although the criti¬ 
cal density depends on the temperature, n cr changes by a small 
percent between 200 K and 500 K. Therefore, we adopt 250 K 
as a high-temperature limit, since the temperature dependence 
of the exponential term in the two-level system approximation 
of N(C + ) is stronger. 

When using a higher temperature and density, the estimated 
column density decreases by less than 20% when considering 
H i or H 2 as the collision partner. The column density of [C 11 ] 
would increase by larger factors if lower densities and tempera¬ 
tures were used. When using electrons as the collision partners, 
N(C + ) saturates at densities above 100 cm -3 , and temperatures 
above 400 K. These results emphasize our point that the val¬ 
ues of N(C + ) obtained with the LTE approximation (or with the 
temperature and density assumed above) should be regarded as 
lower limits, whether the LTE conditions are met or not in all the 
regions mapped. 

The channel maps of the residual [C 11 ] column density (not 
associated with the dense and halo molecular gas traced by C 18 0 
and [C 1 ]), estimated with Eq. ([2]), is shown in Fig. 14 The chan¬ 
nel maps between 19 kms -1 and 22 kms” 1 are the most af¬ 
fected by the subtraction of the [C 1 ] and C 18 0 emission, thus 
confirming the strong spatial association found in the central ve¬ 
locity range, as shown in Sect. |5.1| The small self-absorption 
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Fig. 13. [C ii] column density enhancement for three collision partners, 
H 2 (top), H I (middle), and e~ (bottom), with respect to the column 
density obtained assuming a temperature of 250 K and a density of 
10 4 cm -3 . This column density is depicted by the contour line equal 
unity. 


that we see will affect the neighboring channels, but not the 
general picture. This is in line with the channel maps from 18- 
20 km s -1 that still follow the structure of M17SW. The column 
density of the [C ii] gas not associated with the star-forming ma¬ 
terial (traced by [C i] and C 18 0) ranges between ~10 14 cm -2 and 
~4xl0 17 cm -2 in the whole region mapped, and among all the 
channels in the 0-40 km s -1 velocity range. 


From the residual (and total) [C n] column density channel 
maps we can also compute the corresponding mass of the gas 
contained in the 24" beam area of the channel maps by assuming 
a gas phase carbon abundance of X(C + /H) = 1.2 x 10 -4 ( j Wake-1 


lam & Herbst 2008, their Table 1) and complete ionization of the 


carbon, according to 


M gas = 1.4m H 


N(\C ii]) 

Xq + /H 


Abeam-) 


(3) 


where Abeam is the area (in cm 2 ) covered by the 24" beam, mn is 
the atomic hydrogen mass (in g), and the factor 1.4 accounts for 
helium and a minor fraction of other heavier elements. 

The velocity distribution of the mass (obtained by adding up 
all the pixels from each velocity channel map created from Eq.[3]) 
is shown in Fig. [15] The gas mass per velocity channel (squares) 
was estimated from the original [C n] spectra. The largest frac¬ 
tion of non-associated gas mass (circles) is found at the higher 
(25-33 kms -1 ) velocity channels. 

When integrating the mass in the 0-40 km s -1 velocity range, 
we find a gas mass of ~4.4xl0 3 M© in the entire region mapped. 
This massjs a facto r ~3 lower than the 1.45xl0 4 M 0 found by 
(Stutzki & Giisten 1990) from C 18 0 observations, which trace 
the non-dissociated cloud core mass. We note that the mass 
obtained from the LTE estimated column density of [C i] (cf. 
Fig. |A.3| ) using the same Eq. [3] is about three orders of mag¬ 
nitude lower than the mass traced by [C n]. This is due to the 
lower intensity and less spatial extension of the [C i] emission 
throughout the mapped region and among all the velocity bins, 
compared to that of [C n]. This is in line with the fact that only 
a small fraction of carbon is expected to be in atomic gas traced 
by the [C i] line. 

Considering the mass estimated from the residual [C n] spec¬ 
tra, the gas mass from the non-associated [C n] emission is 
~2.8xl0 3 Mq. Thus, the estimated gas mass not associated with 
the star-forming material traced by [C i] and C 18 0 corresponds 
to ~64% of the total gas mass traced by the original [C n] emis¬ 
sion. This still amounts to at least 19% of the C 18 0 mass reported 
by |Stutzki & Gtiste~n ( 1990| ). 

A source of uncertainty to consider in our analysis is that 
we are assuming that the volumes of gas corresponding to each 
velocity channel are associated. In other words, a spatial associ¬ 
ation of the [C ii] emission with [C i] and C 18 0 emitting regions 
in the plane of the sky, does not ensure that they are really asso¬ 
ciated along the line of sight. However, the probability that the 
spatial association that we see in the central components of the 
1 km s -1 channel maps is just a projection effect is minimal since 
we are using a relatively narrow velocity width. Hence, the frac¬ 
tion of non-associated mass quoted above should be considered 
a lower limit. 

Since there is no evidence of fast shocks in M17 SW, other 
mechanisms that can produce [C n] emission with larger veloc¬ 
ity dispersion than the denser molecular gas must be considered. 
For instance, the interaction with winds and outflows from the 
ionizing stars can lead to substantial excitation of the [C n] emit¬ 
ting gas. Hence, ablation (e .g.,|Castor et al.|1975[|Weaver et al 


1977 |Tenorio-Tagle||1979} |Henley et al.||2012[ ~and references 


therein) and, probably slow shock- interaction due to radiative 
pressure (e.g., |Goodwin||1997[ |Krumholz et al.|[20l0t |Dale & 


|Bonn ell 2011) and references therein), have to be considered to 
model and interpret the [C n] emission not associated with the 
star-forming material in M17 SW. However, [C n] may also be 
present in extended low density gas around H n regions produced 
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Fig. 14. Velocity channel maps at 1 km s 1 width of the residual [C n] column density (cm 2 in logio scale), estimated assuming LTE conditions, 
that is not associated with the dense and halo molecular gas traced by C 18 0 and [C i], respectively. 



Velocity [km/sj 


1-1 <1*S 




Fig. 15. Gas mass (squares) estimated from the [C n] 158 pm emission 
with Eq. { 3 } at each velocity channel in the range 0 km s -1 to 40 km s -1 . 
The gas mass not associated in the spatial distribution of each velocity 
channel map with star-forming material traced by the C 18 0 / = 2 —> 1 
and [C 1 ] 609 pm lines is shown by filled circles. 


by far ultraviolet (FUV) photons (£>11.26 eV for the first ion¬ 
izing potential of atomic carbon) from the ionizing stars. With 
such high energy photons, low H 2 densities also mean higher 
abundance (i.e., densities) of H 1 and free electrons e~, which 
would then become equally important collision partners of [C 11 ], 
with critical densities of about 3xl0 3 cm -3 and 10 cm -3 , respec¬ 
tively, at 250 K (the higher the temperature, the higher the criti¬ 
cal densities). Thus, compensating the lower densities of H 2 and, 
most likely, keeping the LTE assumption for [C 11 ] valid. There- 
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fore, warm C + could be present at lower densities (and higher 
gas temperatures) than assumed. If so, then extreme ultravio¬ 
let (EUV) photons should also produce [N 11 ] emission if their 
energy is larger than 14.5 eV (the first ionizing potential of ni¬ 
trogen) and lower than 24.38 eV (the second ionizing potential 
of carbon). In such a zone of EUV photons, the [C 11 ] and [N 11 ] 
emission should co-exist, and show a high degree of spatial asso¬ 
ciation. This can be checked observationally, using the GREAT 
instrument on board SOFIA. Zones with photon energies larger 
than 24.59 eV (the first ionizing potential of helium) accounts 
for at least 10% of the gas in Ml7 SW (depending on the spatial 
resolution of the observations) as estimated from observations of 


the He + /H ratio (e.g., Peimbert et al. 1988 Tsivilev & Krasnov 
\1999) . Between the two extremes of molecular and ionized hy¬ 
drogen, the [C 11 ] emission can a lso co-exist with neutral atomic 
hydrogen H 1 , as shown by ( [Brogan & Troland|2001] ) in Ml7 SW. 
We discuss these gas phases in more detail in the following sec¬ 
tion. 


5.3. [C 11 ] in the three gas phases 

As discussed above, there are basically three different regimes 
that contribute to the [C 11 ] emission; the highly ionized gas 
where electrons dominate (H 11 ), the atomic hydrogen layer (H 1 ), 
and the molecular hydrogen gas (H 2 ) suffused with sufficient UV 
to keep CO dissociated and to ionize neutral carbon efficientl}]^] 
Following the method described above, we include the high res- 


8 There is also the CO-dark molecular gas |Wolfire et al. 2010) that we 
trace in part through our [C 1 ] line. However our [C 1 ] sensitivity is not 
good enough to also measure the more diffuse molecular, CO-dark gas. 
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Fig. 16. Fraction of the [C n] 158 fim emitting region correlated (at each 


1 km s channel) with the optical depth of HI (from Brogan & Troland 
|2001| ). The corresponding correlation coefficient is overlaid. 
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Fig. 17. Residual [C n] 158 /um spectrum (dashed line and gray filled 
histogram) at offset position (-30",-15"), after subtracting, from the 
original [C n] spectrum, the scaled up (*) spectra of model l-(top), 2- 
(middle) and 3-(bottom), as well as the optical depth of H I, r(H i) (see 
text). All negative (noise) channels in the scaled up and residual spectra 
are set to zero. The residual [C n] spectrum is shifted in -1 K for clarity. 


oluti on VLA map of the veloci ty-resolved optical depth, r(H i), 
from |Brogan & Troland ( |2001[ ), convolved with a 24" beam and 
re-sampled to 1 kms” 1 channel width. The channel-by-channel 
spatial correlation between [C n] emission and r(H i) shown in 


Table 1. Fraction of average [C n] emission in the three gas phases. 


Model No. 

H 11 fraction 

H 1 fraction 

H 2 fraction 

Model l a 

42.7% 

20.9% 

36.4% 

Model 2 b 

36.2% 

16.8% 

47.0% 

Model 3 C 

36.8% 

17.4% 

45.8% 


Notes. (a) Combination: [C n] - {r(H i) + [C i] 609 jum + C 18 0 (2 - 1)} 
(b) Combination: [C n] - {r(H i) + [C i] 609 jjm + 12 CO (1 - 0)} (c) Com¬ 
bination: [C II] - {r(H I) + 12 CO (1 - 0) + C 18 0 (2 - 1)} 


Fig. [IT] indicates that most of the [C n] emission associated with 
the H i gas is found at the lower (< 20 km s -1 ) velocity channels. 

In the previous section we estimated the [C n] column den¬ 
sity and hydrogen mass of the gas not associated with the 
relatively compact and dense star-forming material traced by 
C 18 0 / = 2 —> 1, and [C i] 609 jam. The [C n] emission asso¬ 
ciated with the entire molecular gas phase, however, also com¬ 
prises the diffuse and more extended H 2 gas. Therefore, we now 
use the 12 CO / = 1 —> 0 line, as the canonical tracer of H 2 . As 
12 CO / = 1 —> 0 is optically thick throughout large parts of the 
map, this will provide a lower limit for the [C 11 ] emission from 
the diffuse molecular material. 

In order to analyze the impact of these line tracers in the 
residual [C 11 ] emission (and hence, the [C 11 ] column and as¬ 
sociated gas mass), we tested three different combinations of 
gas tracers: (1) r(H 1 ) + [C 1 ] 609 jam + C 18 0(2-l); (2) r(H 1 ) 
+ [C 1 ] 609 n m + 12 CO(1-0); and (3) r(H 1 ) + 12 CO(1-0) + 
C 18 0(2-l). Model (3) is included for comparison, since C 18 0 
can complement 12 CO in regions where 12 CO is optically thick. 
The residual [C 11 ] spectrum at offset position (-30",-15"), ob¬ 
tained after subtracting the scaled-up spectra of the three com¬ 
binations mentioned above, are shown in Fig. [IT] (from top to 
bottom). When subtracting (channel by channel) the maximum 
of the three synthetic lines from the original [C 11 ] spectra, we 
obtain the residual [C 11 ] emission that is mostly associated with 
the H 11 regime. It can be seen in the three cases that most of 
the residual [C 11 ] emission is contained in the higher velocity 
channels. Using only the molecular gas tracers (i.e., excluding 
r(H 1 )), we can obtain a second residual [C 11 ] spectra that would 
contain the [C 11 ] emission associated mostly with the H 11 and 
H 1 gas. From these two residual spectra we can then estimate 
the residual [C 11 ] emissions, and their respective column densi¬ 
ties using Eq. ([2]), associated with the three gas phases, according 
to the procedure described in Appendix [E] 

In Table [T] we summarize the fraction of the average over 
the region mapped [C 11 ] emission associated with the three gas 
phases as obtained from the three combinations of line tracers. 
Using [C 1 ] 609 /im and 12 CO / = 1 —> 0, combined with r(H 1 ), 
yields practically the same result given the uncertainties as when 
using C 18 0/ = 2—>1 instead of [C 1 ] 609 /im. This results 
from the high correlation observed between [C 1 ] and C 18 0, al¬ 
though it is not a 1:1 match (probably sensitivity limit driven) as 
shown in Fig. [4] Thus, the slightly higher (~1%) fraction of [C 11 ] 
emission associated with H 2 gas found with model (2), compared 
with that of model (3), may indicate that [C 1 ] traces at least part 
of the CO-dark molecular gas ( Wol fire et al.||2010| ). Therefore, 
we chose model (2) as the most complete one, tracing all the gas 
regimes where [C 11 ] emission can be found. 

The velocity distribution of the fraction of averaged [C 11 ] 
emission associated with the three gas phases obtained with 
model (2) is shown in Fig. [18] The [C 11 ] emission associ- 
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Fig. 18. Fraction of the average over the region mapped residual [C n] emission corresponding to the three gas phases: H n (left), HI (middle), and 
H 2 (right) as obtained with model (2) in Table [T] 


ated with Hi gas is mostly contained in the lower velocity 
(<20 kms -1 ) channels, while the [C n] emission associated with 
the ionized H n gas is contained mainly at the higher veloc¬ 
ity bins (>25 kms -1 ), although part of it is also found in the 
<20 kms -1 velocity range. The central velocity channels (15- 
30 km s -1 ) contain most of the [C n] emission associated with the 
molecular H 2 gas. The corresponding velocity-channel maps of 
the [C 11 ] emission associated with H 11 , H 1 , and H 2 , are shown in 
Fig. |E.1| These channel maps show that the [C 11 ] emission, and 
therefore column density, associated with the ionized gas, peaks 
at the northeast comer of the mapped region, which coincides 
with the position of the ionizing sources (cf. Fig [5]). The fraction 
of [C 11 ] column density (or [C 11 ] emission) associated with the 
molecular gas regime is about 11% larger than the fraction (36%) 
found in the dense star-forming material. This is expected since 
the 12 CO emission has a broader line profile and is also spatially 
more extended than C 18 0 and [C 1 ] (cf. Fig. [ 2 ] and Fig. [Tl]). 

We note that this method has uncertainties, and it gives 
only a first order approximation of the [C 11 ] emission associ¬ 
ated with the three different regimes. The results presented in 
Table [T] should not be taken as a sharp distinction between the 
three gas regimes, since in reality the three gas phases can be 
mixed throughout the region mapped (we elaborate on this in 
the next sections). In particular, the [C 11 ] emission associated 
with the atomic H 1 gas has a large uncertainty as the optical 
depth r(H 1 ) is saturated over a significant part of the region 
along the molecular ridge. The saturated values were replaced 
by a lower limit, according to the continuum and rms level of the 
VFA spectra, as described in detail by Brogan & Troland ( 2001, 
their Sect. 3.3.2). Furthermore, r(H 1 ) obtained from H 1 in ab¬ 
sorption traces atomic hydrogen in the foreground relative to the 
free-free emission from the H 11 region. This might introduce a 
bias since r(H 1 ) traces only the cold H 1 gas, while part of the 
warmer atomic hydrogen can be mixed with the H 2 and H 11 gas 
phases. However, we consider that the warm mixed H 1 gas can 
be at least partially accounted for by the [C 1 ] and 12 CO lines. 
Nevertheless, this is another uncertainty in our method. 

The gas masses associated with the three gas phases could 
be estimated with Eq. 0, by using abundances (and the corre¬ 
sponding mass of atomic and molecular hydrogen) of ionized 
carbon relative to the dominant hydrogen phase, i.e., X(C + /H°), 
X(C + /H 2 ). However, these abundances are not really known, and 
although they could be estimated from a clumpy PDR model, 
which would be the best model currently available for M17 SW 
because of its highly clumpy structure, the uncertainties in the 
values obtained for each gas phase would be very high be¬ 
cause the relative abundances strongly depend on the number 
of clumps, clump sizes, and ambient conditions of each clump, 
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which we have not yet been able to constrain for M17 SW. In ad¬ 
dition, knowing the actual density and temperature of the domi¬ 
nant collision partners in each gas regime would allow us to es¬ 
timate (with non-FTE radiative transfer models) the actual [C 11 ] 
emission associated with the three gas phases. We expect to es¬ 
timate all these parameters in a follow up work. 


5.4. Comparison with radio recombination lines 


Traditionally radio recombination lines (RRF) are considered 
probes of gas conditions within ionized gas. However, the princi¬ 
pal quantum levels are generally not in thermal equilibrium, and 
there are different competing mechanisms governing the line in¬ 
tensities (e.g., maser amplification of background radiation and 
weakening due to underpopulation of the upper quantum lev¬ 
els relative to FTE conditions) as well as the broadening of the 
line widths (e.g., thermal, turbulence and Stark effects - Griem 
|1967| . Furthermore, the line shapes are the result of emission 
from regions with different ambient conditions along the line of 
sight, since the emitting gas is usually optically thin in the cen¬ 
timeter and millimeter regimes. All these effects make of RRFs, 
in the centimeter wavelength range, more ambiguous probes of 
the ambient conditions in ionized gas than originally thought. At 
millimeter wavelengths, however, the Stark broadening is negli- 
gible (AV ~ c(-mrh^) 5/3 <0.1 kms -1 for v > 90 GHz, cf. [Gordon 
|& So rochenko 2002), leaving only thermal and turbulent broad- 
ening at play, and the RRFs tend to be optically thin and do not 
suffer significant departures from FTE conditions (for a review 
see, e.g., |Gordon|1988[[Gordon & Sorochenko|2002 ). 

In our broadband (32 GHz IF bandwidth) OTF maps ob¬ 
tained with the IRAM 30m telescope (cf. Sect. [2]), we detect a 
number of hydrogen RRFs in the 3mm band: H39 a, H40 a, and 
U41a lines at 106.74 GHz, 99.02 GHz, and 92.03 GHz, respec¬ 
tively (the beam size of the H41 a line is 2873). All of them show 
very similar spatial distribution and line shapes. We improved 
the S/N by convolving the maps with a 30" beam. Because of 
our short integration times no other recombination lines were 
detected, preventing further analysis of the excitation conditions 
of the ionized gas at our high spatial resolution. Our goal now is 
simply to verify if the spatial distribution of the [C 11 ] emission 
can be associated with the ionized gas traced by the hydrogen 
recombination lines. 


Figure 19 shows the H39 a, [C 11 ], and 12 CO / = 1 —> 0 
integrated intensities, as well as the r(H 1 ) optical depth over¬ 
lay (right panel) on the H41 a map. The hydrogen recombina¬ 
tion lines follow a relatively homogeneous distribution about 45° 
along the ionization front, and they peak at the M17-UC1 ultra¬ 
compact region, which is also surrounded by a number of em- 
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Fig. 19. Left - Velocity-integrated intensity maps of H41cr (gray), H39cr (red contour), [C n] 158 /dm (green contour), and 12 CO / = 1 —» 0 (blue 
contour). The contour lines (from t hin t o thick) are 10% (d ashed line), 25%, 50%, 75%, and 90% of the respective peak emissions. The stars 
indicate the O and B ionizing stars (Beetz et al. 1976, Ha nson et al.|1997|>. The refe rence position (A a = 0, AS = 0), marked with a cross, is as in 
Fig. [2] The ultracompact H 11 region M17-UC1 and four H 2 Q masers ( [Johnson et al.|1998| ) are marked by the circle and “+” symbols, respectively. 
The small purple circles correspond to the heavily obscured (E median > 2.5 keV, A v > 10 mag) population of X-ray sources around the M17-UC1 
region (Fig. 10 in Broos et al. 2007; coordinates from the VizieR catalog). Right - Same as on the left, but with r(H 1 ) instead of H39 a. All maps 
have been convolved with a 30" beam, to increase the S/N of the H41of map. 


bedded X-ray sources (Broos et ah] 20071 ). We note that there 
are many more X-ray sources reported by Broos et al., several 
of them with a stellar counterpart, but for clarity of the fig¬ 
ures we show here only those associated with the M17-UC1 
region. The H 1 optical depth follows the distribution of the 
RRLs, but their peak emissions are not correlated. In fact, the 
velocity integrated r(H 1 ) peaks around the same position as the 
12 CO / = 1 —> 0 line intensity. The [C 11 ] emission, instead, cov¬ 
ers the entire region mapped, and its intensity peaks between the 
ionized gas traced by H41 a and the molecular gas traced by the 
12 CO / = 1 -> 0 line. 

Figures |20| and [2T| show the spectra of the H41 a, [C n], and 
12 CO J ~ 0 lines at two offset positions close to the ion¬ 

izing sources, (30",100") and (20",70"), and along a strip line 
at P.A.=63°. Although fainter at the northeast region (Fig. [20|, 
the H41 a line shows emission only at the higher velocity range 
(>20 kms -1 ) coinciding with the velocity channels of the [C 11 ] 
lines found to be associated mostly with the H 11 gas in Sect. |5.3| 
At the southwest region (Fig. \2A\ the H41 a line is stronger, but 
its shape is asymmetric and much broader than any of the other 
lines we observed in Ml7 SW. It is hard to associate its line 
shape to any of the line structures observed in the [C 11 ] line nor 
in r(H 1 ). Since the H41 a line is expected to be optically thin 
(i.e., no optical depth nor self-absorption effects), its line shape 
is most likely formed by several layers of ionized gas (with dif¬ 
ferent Visr) along the line of sight, with each velocity compo¬ 
nent being affected by a combination of thermal and turbulence 
broadening. 

From Fig. [18] we know that the velocity ranges of the 
peak [C 11 ] (residual) emission found to be associated with the 
Hu, Hi, and H 2 gas, are 32-33 kms -1 , 9-10 kms -1 , and 21- 
22 kms -1 , respectively. Figure 22 shows the corresponding 
channel maps of H41or, the original [C 11 ] emission, r(H 1 ), and 
the 12 CO / = 1 —> 0 line. We note that the spatial distribu¬ 
tion of the H41of line does not change significantly over these 
three different velocity ranges. At 32-33 kms -1 , the [C 11 ] emis¬ 
sion peaks at the northeast region, around the ionizing sources, 
while the H41 a emission still peaks at the M17-UC1 region (cf. 


velocity range 9-10 kms -1 , identified mostly with the H 1 gas 
regime, the emission from all the lines seems to emerge mostly 
from the region of the ionization front traced by the H41 a line, 
indicating that the three gas regimes must be mixed at these ve¬ 
locities. The fact that the ionization and molecular dissociation 
front are mixed, can be explained if advection is considered (e.g., 
Bertoldi & Draine||1996| ). This dynamical process would make 
the ionized region larger compared to the atomic region, and ex¬ 
tended towards the molecular region. This, in turn, will lead to 
a larger contribution of the ionized gas to PDR diagnostics like 
[C 11 ], as shown in models by |Abel et al.[ ( [2005] ). 

As previously noted by Bro gan et al.| ( |1999| ), and from 
Figs.[2Q|and[2H the velocity structures of H 1 and [C 11 ] are differ¬ 
ent. A narrow component of shocked H 1 gas (streaming towards 
us) in M17 was suggested by |Brogan et al.| ( p~999| , in the range 
11-17 kms -1 , where we also find [C 11 ] and '“CO components, 
but only at the edge of the ionization front (cf. Fig. [21] offset po¬ 
sition (-30",-15")). In the 9-10 kms -1 velocity range, which 
is associated mostly with atomic gas, the H 1 optical depth peaks 
along the H41 a line since r(H 1 ) is estimated from 1 in absorp¬ 
tion. Actually, this channel map is part of the prominent r(H 1 ) 
velocity component in the 0-11 kms -1 range (cf. Fig. 20), that 
has corresponding associated [C 11 ] and 12 CO / = 1 —> 0 emis¬ 
sion. This component in r(H 1 ) was not mentioned in the original 
papers by |Brogan et al. ( |1999| ); [Brogan & Troland| ( [2001| ). This 
can also be shocked ionized ([C 11 ]), atomic (H 1 ), and molecular 
(CO) gas streaming towards us. 

The strong JG 11 ] velocity component found around Vi sr = 


32 km s (Fig. 20 top two panels) has no evident counterpart in 


Fig. 19). At 21-22 km s , the H41 a peak is at the southeast re¬ 


gion of the ionization front, while the distribution of the [C 11 ] 
emission follows closely that of the 12 CO / = 1 —> 0. In the 


the other atomic or molecular tracers. These differences, in par¬ 
ticular with H 1 , may arise because the [C 11 ] 158 jam emission 
tends to trace warmer gas than the H 1 absorption, which traces 
cold gas. The faint H41 a emission around Vi sr = 32 kms -1 , 
barely detected with S/N-3 in our OTF maps at offset position 
(20",70"), shows a line width (FWHM) of AV = 20 ± 3 km s -1 
from a Gaussian fit. This line width is consistent with the thermal 
broadening (A VG-t hermai ~ cx 7.16233_x 10 -7 (-^) 1/2 , with T in K 
and M in amu, cf. Gordon & Sorochenko|2002[ their Eq. 2.22) 
expected for the LTE elect ron temperature T * = 10700 ± 700 K 
estimated towards M17 by |Gordon| ( 1989| ). At offset (30",100") 
we find a S/N<3 for the H41 a line. At these positions, (20",70") 
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Fig. 20. Spectra of several lines observed at approximated (±3") offset 
positions along the strip line at P.A. 63° (A£ = Aa/2) and at positions 
close to the ionizing stars. All the spectra have been resampled to a 
1 km s- 1 resolution and convolved with a 30" beam, to increase the S/N 
of the H41 a map. A Gaussian was fit to the H41cr line at the top two 
positions. 


and (30",100"), the [C n] component around Vi sr = 32 kms' 1 
have a Gaussian fit FWHM of ~ 5.5 ± 0.5 km s -1 , also consistent 
with the thermal broadening expected for the T* quoted above. 
This indicates that the [C n] emission around Vi sr = 32 km s -1 is 
associated with warm gas, and probably shocked by the proxim¬ 
ity of the ionizing sources, streaming away from us and from the 
bulk of the molecular unshocked gas found around the 20 km s -1 
component. Higher S/N maps of the H41 a and other RRLs are 
needed to confirm the detection of ionized gas at these posi¬ 
tions. Maps of [N n] and the hydrogen recombination lines (in 
the >3 mm wavelength range), tracing the transitions between 
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Fig. 21. Spectra of several lines observed at approximated (±3") offset 
positions along the strip line at RA. 63° (A5 = Aa/2). All the spectra 
have been resampled to a 1 km s -1 resolution and convolved with a 30" 
beam, to increase the S/N of the H41 a map. 


medium and large principal quantum numbers, and even the 
fainter carbon recombination lines, may help explain the nature 
and ambient conditions of the 32 km s _I component found in the 
[C n] emission. 


5.5. Comparison with other sources and implications for 
extragalactic studies 

A study of the H n region S125 showed that up to 40% of the 
[C n] 158 pm line in tensity, observed with ISQ -SWL, arises 
from the H n region ( [Aannestad & Emery) ||2003| . Comparing 
and modeling the [C n] and [N n] emission obtained from Her- 
schel/PACS, |Bernard-Salas et~aL] ( |2012| ) found that at most 18% 
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Fig. 22. Channel maps of the H41 a line (gray background, in Kkm s _1 ) 
with overlays of the [C n] (red contours), the r(H i) optical depth (blue 
contours) (left panels ), and the 12 CO / = 1 —> 0 (green contours) (right 
panels ), at the velocity ranges of the peak [C n] emission associated 
with (from top to bottom) the H II, H I, and H 2 gas regimes. The corre¬ 
sponding velocity ranges are shown in the top-left of each channel map. 
The contour lines (from thin to thick) are 10% (dashed line), 25%, 50%, 
75%, and 90% of the respective peak emissions. All the maps were con¬ 
volved with a 30" beam to increase the S/N of the H41cr map. 


of the [C n] emission in the Orion Bar originated in the H n re¬ 
gion. This is a factor of -2 lower than found for SI25, and a 
factor of - 3.6 times lower than what we estimate for M17 SW. 
However, this is not a one-to-one comparison, since their study 
is focused in a relatively small scale, covering mostly the dense 
PDR in the Orion Bar. Hence, most of the actual H n region 
(where more [C ii] emission than quoted by them can be present) 
is not seen in their data set. Therefore, large scale maps of [C n] 
and other ionized, neutral, and molecular tracers like the ones 
we present are relevant not only in order to understand the full 
picture of the H n-PDR boundary, which is more complex than 
predicted by the classical 1-D models, but also to understand 
and properly interpret the observations of these diagnostic lines 
in extragalactic sources. 

As m entioned at the end of Sect. |5.2| and as shown in 
Sect. |5.4| the [C n] line can originate both in the PDR (molec¬ 
ular and neutral atomic gas) and in the H n region. When ob¬ 
served toward other galaxies, the [C n] emission collected in 
one beam (typically covering large scales), will naturally come 
mainly from several unresolved giant molecular clouds and H n 
regions, with different orientations with respect to the imping¬ 
ing radiation field along the line of sight. Thus, it is important 
to characterize in detail the [C n] contribution from these differ¬ 
ent environments. The [C n] emission is important in extragalac¬ 
tic studies for redshift determinations, and also to estimate the 
star formation rate (SFR) from the [C n] luminosity (e.g., |Stacey| 
et al.|199T||Meijerink et al.|2007HLuhman et al.|2003| ). 

Recent results from 130 galaxies observed with Her- 
schel/PACS by |Sargsyan et al.| ( |2014| ) confirm that [C n] traces 
the same starburst component of sources as measured with mid- 
infrared PAH, neon emission line diagnostics, and bolometric 
luminosities of reradiating dust from starbursts. By using a mod¬ 
ified Cloudy code [Abel et aL| ( j2005| ) estimated that about 30% 
of the [C n] emission in the starburst galaxy NGC 253, ob¬ 
served with the KAO ( [Carral et al.|1994| , arise from the ionized 


medium. Using PACS data to map the FIR emission of the [C n] 
line in the spiral galaxy M3 3, and models of photo-ionization 
and photon-dominated regions, |Mookerjea et al.| ( |2011| ) found 
that between 20% and 30% of this emission comes from the H n 
region. 

Compared with our lower limit estimates, these results for 
extragalactic sources are at least a factor of two lower than the 
values we find for M17 SW when considering only the [C n] 
emis sion not associated with the dense star-forming material (cf. 
Sect. |5.2| ). They are more than 6% lower than what we estimate 
for the [C n] emission associated mostly with the H n region (cf. 
Sect. H3>. We believe that there are three main reasons for these 
differences: 


(1) Their models of static geometries do not consider that their 
beams collect emission from many sources with diverse ge¬ 
ometries (from edge-on to face-on, and many other orien¬ 
tations in between), systemic velocities, highly structured 
medium (that would allow FUV radiation to permeate the re¬ 
gion ionizing and heating the gas on larger spatial scales), 
and dynamical processes like advection (as mentioned in 
Sect. |5.4| ) that would increase the contribution from the ion¬ 
ized medium. 

(2) The models assume that the gas is in pressure equilibrium, 
so the equation of state is dominated by gas pressure since 
turbulent and magnetic pressures are not included. If turbu¬ 
lent or magnetic pressure dominate instead, as is the case 
of M17 SW (Pell egrini et al.||2Q07] ), then the density law 
and hence the collisional excitation of the PDR diagnostic 
lines throughout the H n and photon-dominated regions will 
be different. 


(3) Their studies are based on the total velocity-integrated 
intensities of the [C n] and other line emissions, pro¬ 
vided by the spectral resolution of the instruments used. 
Whereas the velocity-resolved [C n] spectra obtained with 
SOFIA/GREAT towards M17 SW, shows that -80% of the 
velocity range covered by the [C n] line is associated with a 
mixture of ionized medium, neutral gas, and diffuse molec¬ 
ular gas, and only -20% of the velocity range is associ¬ 
ated mostly w ith t he den se (star-forming) molecular mate¬ 
rial (cf. Sects. [5H] to |5.4| ). While this result is particular for 
Ml7 SW, we have found similar results in other regions (e.g., 
NGC 3603, Perez-Beaupuits et al][, in prep.]), showing that 
this characteristic should be taken into account in PDR mod¬ 
eling of Galactic and extragalactic sources. 


All the reasons mentioned above lead us to think that the 
fraction of [C n] emission coming from the H n regions have 
been underestimated in previous studies. Their conclusions re¬ 
garding the fraction of [C n] emission emerging from ongo¬ 
ing star formation, and their estimates of SFRs in comparison 
with other SFR diagnostics, may change. In order to judge how 
much more valid our estimate is, models of photo-ionization and 
photon-dominated regions that include dynamical processes and 
magnetic pressure are needed. In addition, a significant number 
of large scale maps of [C n]and other diagnostics of ionized gas 
toward Galactic sources are required to increase the statistical 
significance of our results. Large scale maps of many sources 
can be achieved with the upGREAT receiver arrays, the second 
generation receivers for the GREAT project, that will make map¬ 
ping the large scale regions more efficient. The upGREAT low 
frequency array (1.9-2.5 THz with 14 pixels) is currently sched¬ 
uled for commissioning onboard SOFIA in May 2015. The up¬ 
GREAT high frequency array (4.7 THz with 7 pixels) will follow 


Article number, page 19 of [27] 




































































A&A proofs: manuscript no. 25020_printer 


about one year later. We expect to show the results of this long 
term study in follow up work. 

6. Conclusions 

We used the dual channel DSB receiver FLASH on the APEX 
telescope to map (with 1277 resolution) a region of about 4.1 pc 
x 4.7 pc in the 3 Pi —> 3 Po 609 pm (J = 1 —> 0) fine-structure 
transition of [C i], towards the star-forming region M17 SW. We 
also used the broadband EMIR receivers on the IRAM 30m tele¬ 
scope to map a similar area of 360"x300" in the 3mm, 2mm, 
and 1mm bands. 

We combine these data with the previously observed and 
published [C i] 3 P^ —> 3 Pi 370 pm and [C n] 158 pm data to 
discuss the physical properties of the dense interstellar medium 
in the M17 SW ridge. 

Because of the complex structure of the M17 SW and the 
availability of velocity resolved spectra, we have performed all 
our analysis on 1 kms -1 -wide velocity channel maps. 

6.1. Excitation and column densities of [C i] 

Combining our earlier observation of the [C i] 3 P 2 —> 
3 Pi 370 pm (J = 2 —> 1) fine-structure line together with 
the new [C i] 3 Pi —» 3 Po 609 pm data, we found that the 
R = 7([C i] 369//m)//([C i] 609/mi) ratio is larger than unity in 
most of the regions mapped, at the central (10-24 kms -1 ) ve¬ 
locity channels where the bulk (>20%) of the [C i] emission is 
found. We estimate the excitation temperature T ex and column 
density of [C i] using an optically thin approximation and a non- 
LTE method. We found that T ex ranges between -40 K and -100 
K in the inner region (i.e., southwest from the ionization front). 
While the [C i] column density ranges between ~4xl0 13 cm -2 
and -10 17 cm -2 throughout the whole region mapped and among 
all the velocity channels. In the region where the [C i] emission 
is > 50% of its peak integrated intensity, column densities up to 
10 17 cm -2 are only reached in the central (17-22 km s -1 ) velocity 
channels. 

6.2. Comparison of [C n], [C i], and CO 

We used our recent SOFIA/GREAT velocity-resolved [C n] map 
to analyze the spatial and velocity association of the [C n] emis¬ 
sion with the diffuse and molecular gas in M17 SW. For that we 
used the [C i] lines from APEX/FLASH and the 12 CO and iso¬ 
tope lines from IRAM30m/EMIR. 

The [C n] emission was found to be associated with the 
other gas tracers in 20%-80% of the mapped region, but only at 
the central (15 kms -1 and 23 kms -1 ) velocity channels, which 
means that only -20% of the velocity range (-40 kms -1 ) that 
the [C n] line spans in our velocity-resolved spectra is associated 
with the star-forming material in M17 SW. 

For the non-associated [C n] gas, in the 1 kms -1 wide 
channel maps we estimated column densities ranging from 
-10 14 cm -2 km -1 s to ~5xl0 17 cm -2 km -1 s across the region 
mapped. The largest fraction of non-associated atomic mass 
contained in the mapped region, is found at the higher (25- 
33 km s -1 ) velocity channels. The total non-associated gas mass 
(integrated over the 0-40 kms -1 range) is ~2.8xl0 3 M 0 . This 
corresponds to a very large fraction (-64%) of the total mass 
(~4.4xl0 3 Mq) traced by the [C n] emission. In other words, 
most of the gas traced by the [C n] emission in M17 SW is not 
associated with the star-forming material. The good match of the 
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[C i] emission with the C 18 0 / = 2 —> 1 map indicates that both 
tracers are seen from the same, star-forming material. 

6.3. Comparison of [C n], H i, and hydrogen recombination 
lines 

When using the optical depth of H i in combination with the 
[C i] 609 pm and 12 CO / = 1 —> 0 maps, we found that the 
[C n] emission associated with H i is contained in the lower 
(<20 km s -1 ) velocity channels, while the [C n] emission associ¬ 
ated with the H n gas phase is found mostly at the higher veloc¬ 
ity (>25 km s -1 ) channels. Most of the [C n] emission associated 
with the diffuse and dense molecular H 2 gas, is found in the cen¬ 
tral velocity channels (15-30 km s -1 ). From our preferred model 
(2), we found that 36.2%, 16.8%, and 47.0% of the [C 11 ] emis¬ 
sion is contained in the H 11 , H 1 , and H 2 regimes, respectively. 

Overlays between [Cn], r(H 1 ), 12 CO / = 1 —> 0, and 
the H41 a line at velocity ranges associated with the three gas 
regimes indicate that the H 11 region is mixed with the atomic 
and part of the molecular dissociation regions, in agreement with 
the highly clumped structure and dynamical processes at play in 
M17 SW. 

Considering M17 SW as a proxy for active galaxies, these re¬ 
sults are also relevant to extra-galactic studies in which [C 11 ] is 
often used as a tracer of star formation rates. We have estimated 
a fraction of the [C 11 ] emission not associated with the dense 
star-forming material, based on large maps of velocity-resolved 
spectra. Our estimates are up to two times larger than estimates 
done in extragalactic sources, based on velocity-integrated [C 11 ] 
intensities. Other Galactic molecular clouds, for which velocity- 
resolved [C 11 ] maps are also available, will be analyzed in a 
follow-up work. With the enhanced observing capabilities of 
SOFIA’S second generation heterodyne arrays, more sources can 
be observed in the future, increasing the statistics. 
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Fig. A.l. Velocity channel maps at 1 km s 1 width of the optical depth 
ri_>o of the [C I] 1 —> 0 line, estimated assuming LTE conditions in 
M17 SW. 



Fig. A.2. Velocity channel maps at 1 km s 1 width of the optical depth 
i" 2 —»i of the [C i] 2 —» 1 line, estimated assuming LTE conditions in 
M17 SW. 
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Fig. A. 3. Velocity channel maps at 1 km s -1 width of the column density 
(cm -2 ) of [C i] (in logio scale), estimated using the excitation tempera¬ 
ture from Fig. [6] 

Appendix A: LTE Analysis of [C i] 

Following |Frerking et al.| ( [1989| and Schneider et al. (2003, their 
Appendix A), the optical depths of the [C i] J = 1 —> 0 and 
J = 2 —> 1 lines can be estimated from the excitation tempera¬ 
ture and the peak intensity of both lines, assuming a beam filling 
factor of unity. Knowing the excitation temperature and the op¬ 
tical depth of the / = 1 —> 0 line, the column density N([C i]) 
can be computed as well. Figures [ATT| and [A2| shows the chan¬ 
nel maps corre sponding to the optical depths of both [C i] lines, 
while Fig. |A.3| shows the channel maps of the [C i] column den¬ 
sity, as estimated assuming LTE conditions. 


Appendix B: Non-LTE Excitation Analysis of [C i] 

Following the procedure described in (Perez-Beaupuits et al. 
|2007( |2009| ), we use the radiative transfer code RADEX ( |van 


der Tak et al .|2007| ) to create a data cube containing the intensi¬ 
ties, as well as the excitation temperatures and optical depths 
of the two [C i] transitions, in function of a range of kinetic 
temperatures 7V, number densities ?z(H 2 ) (i.e., excitation condi¬ 
tions), and column densities per line width N/AV. The collision 
rates used were taken from the LAMDA database (Schoier et al. 
12055). 

We only used collisions with H 2 since this is the most abun¬ 
dant molecule. Other collision partners can be H and He. Al¬ 
though their collision cross sections are comparable, H 2 is about 
5 times more abundant than He, and H is at least one order of 
magnitude less abundant than H 2 in the dense cores of molecu¬ 
lar clouds (e.g., |Meijerink & Spaans|2005] ). Hence, including H 
and He as additional collision partners would not produce a sig¬ 
nificant change in our results. We also assumed a homogeneous 
spherical symmetry in the clumps for the escape-probability for¬ 
malism. 

The original RADEX code was modified to include dust 
background emission as a diluted blackbody radiation field, as 
in |Poelman & Spaans| ( |2005| ) and |Perez-Beaupuits et al.| ( [2O09 ). 
The total background radiation is modeled as a composite be¬ 
tween the cosmic background radiation (CMB), assumed to be a 
blackbody function at 2.73 K, and the diluted dust radiation es¬ 
timated as Tdust x B(Jdust\ where B(Td us t ) is the Planck function 
and the dust continuum optical depth TdustW is defined by |Hol- 


lenbach et al. (1991) as TdustW - Doo//m(100/rni//l). We adopted 
an average dust background temperature Tdust = 50 K and the 


high FIR opacity Tiooyum = 0.106 found by|Meixner et al.|(|l992|) 
in M17SW. 

We assume that the emission collected by the beam has a 
homogeneous elliptical Gaussian distribution and that the cou¬ 
pling factor of our beam to the source distribution is unity, so we 
can compare directly with the output of RADEX, which is the 
Rayleigh-Jeans equivalent radiation temperature T R emitted by 
the source. 

We explored all the possible excitation conditions within the 
given range that can lead to the observed radiation temperatures 
and the line ratios between the two [C 1 ] lines. The line ratios 
and the peak temperature of the lower-/ line involved in each 
ratio were used to constrain the excitation conditions. Including 
the rms of the observed spectra and uncertainties in all the as¬ 
sumptions mentioned above, a 20% error of the ratios and peak 
temperatures are used to define a range of values for TV, ft(H 2 ), 
and N/AV within which the RADEX output is selected as a valid 
solution. 

The volume density explored ranges between 10 2 cm -3 and 
10 4 cm -3 , the kinetic temperature varies from 10 K to 500 
K, and the column density per line width lies between 10 10 
cm -2 km -1 s and IO 20 cm -2 km -1 s. Since the optical depth, and 


the line intensities, are proportional to the column density per 
line width N/AV , we generated the RADEX data cube assuming 
the AV = 1 kms' 1 of the velocity channels in the re-sampled 
spectra. In order to constrain the solutions, we fit the line ra¬ 
tio between the peak temperatures of the transitions and the 
radiation temperature of the lower transition (/ = 1 —> 0), 
which are R ~ 1.08 and T R « T mb = 31.9 K at offset posi¬ 
tion (-130",-10"), and R ~ 1.06 and TV ~ T mb = 22.4 K at 
offset position (-130",-70"). 

Figures B to |B.2| show gray scale and contour maps 
of the excitation conditions found to reproduce the observed 
[C 1 ] line ratios and peak temperatures at the offset posi¬ 
tions (-130",-10") and (-130",-70"), respectively. The val¬ 
ues shown correspond to the average of all the possible N/AV , 
T ex , and r found for each pair of excitation conditions (ft(H 2 ) 
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Fig. B.l. Excitation map (top panel ) for the [C I] line ratio observed 

at offset position (-130",-10") in the velocity channel 19-20 kms -1 . 
The contours and labels correspond to the column density per line width 
N/AV (cm -2 km -1 s, in logio scale). In the middle panel the contours 
and labels correspond to the excitation temperature T ex (K) of the J - 
1 —> 0 (left) and / = 2 —> 1 {right) transitions, while the bottom panel 
shows the corresponding optical depths for each line. 

and T k ). The curvature in the solutions depict the dichotomy be¬ 
tween the kinetic temperature and the density of the collision 
partner. In other words, solutions for the observed values can be 
found for higher temperatures and lower densities, but also for 
lower T k and higher n{ H 2 ). The column density per line width 
does not change significantly along the solution curve, but it 
does change slightly across the curves, and especially at lower 
(< 100 K) kinetic temperatures. This means that for a given ki¬ 
netic temperature, N/AV will show a small variation in function 
of density. 

Appendix C: Correlation between line tracers 

In order to verify the apparent spatial correlation between two 
species, we use the correlation between a specific pixel value of 
a map, and the same pixel in the map of another tracer. This can 
be done only in the maps with the same dimensions and spa¬ 
tial resolutions. Therefore, we first convolved all the maps to 
the larger beam size (24") of the C 18 0/= 1-0 line, and we 
used the SOFIA/GREAT map of the [C 11 ] 158 pm line (which 
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Fig. B.2. Excitation map {top panel) for the [C I] line ratio observed 
at offset position (-130",-70") in the velocity channel 19-20 kms -1 . 
The contours and labels correspond to the column density per line width 
N/AV (cm -2 km -1 s, in log i 0 scale). In the middle panel the contours 
and labels correspond to the excitation temperature T ex (K) of the J = 
1 —» 0 {left) and / = 2 —> 1 {right) transitions, while the bottom panel 
shows the corresponding optical depths for each line. 

covers the smaller region) as a template to create all the [C i] 
and CO maps. The sample correlation coefficient commonly 
used to estimate r ^ between the images X and Y is the Pear¬ 
son’s product-moment correlation coefficient defined by ( |Pear-| 
son |1920[ |Rodgers & Nicewander|1988| ), 

_ cov[X,Y} _ ZiZjlXij-XiiYjj-Y] 
y/vc,r\X\var\Y\ Yffij ~ X] 2 Z; ffjij ~ ?] 2 ' 

(C.l) 

where X t j and Y t j are the pixel values of two given maps or 
images (e.g., [C n] and 12 CO / = 1 —> 0), and X and Y are 
the average pixel value of the respective maps. We note that 
the uncertainty of the correlation coefficient can be approxi¬ 
mated as <r r ~ (1 - rl y )/ yjn - 2, and since we have maps with 
n = 31x18 = 558 pixels, the uncertainty will always be between 
10 -2 and 10 -3 . This correlation is applied to the velocity inte¬ 
grated intensity maps, and to the 1 kms -1 width channel maps. 
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We note that because of the positive correlation between the in¬ 
tensities of different tracers, r xy ranges between 0 and 1. As proof 
of concept, the scatter plot and the associated correlation coef¬ 
ficient between the ve locity -integrated intensity of several line 
tracers is shown in Fig. |C.l The application of the scatter plot to 
the channel maps is shown in Fig. |C.2| for two test cases. We note 
that the correlation between the velocity-integrated intensity, as 
well as the channel maps, of the / = 1 —> 0 and / = 2 —> 1 tran¬ 
sitions of 12 CO is very good, even at the faintest emission of the 
lower velocity (> 3.5 kms -1 ) channels, but the good correlation 
is lost at the higher velocity (> 28.5 km s -1 ) channels. In the case 
of the [C i] 370 pm (X-axis) versus 13 CO J = 2 —> 1 (Y-axis), the 
pixel values are more scattered at the lower (< 16 kms -1 ) and 
higher (> 22 kms -1 ) velocity channels because the line intensi¬ 
ties are fainter and, hence, those channel maps are more affected 
by noise. 


Appendix D: Spatial association in channel maps 

We use a simple method to estimate the fraction of the re¬ 
gion mapped where two emission lines are associated in nar¬ 
row (1 kms -1 ) width channel maps. All maps are first con¬ 
volved to the lowest spatial resolution (24" FWHM beam) of 
the C 18 0 J = 1 —>0 map in order to increase the S/N, and the 
size of the [C i], 12 CO and isotope maps is limited to the region 
mapped in [C n] by using the [C n] data cube as a template in 
GILDAS/CLASS. In this way, we produce spectral cubes with 
the same dimensions and number of pixels. 

Then we determine which region of the maps have signif¬ 
icant emission. This can be done naturally by using the rms 
(noise) level of the spectra corresponding to each pixel in the 
map (i.e., a 3-cr level detection), or by defining a threshold for 
the intensities. Since our maps (especially that of the [C n] line) 
are not homogeneously sampled, the rms level varies among the 
different pixels. Hence, we prefer to use a fraction of the global 
peak (maximum) integrated intensity found among all the chan¬ 
nel maps as a threshold. This provides a unique value that is 
used to determine whether the emission of some region (pixel) 
in a particular channel map is significant or not. 

We create binary images for each channel map by assign¬ 
ing a zero to all the pixels with intensity values lower than the 
threshold, and a value of unity to all the pixels that have intensity 
values larger than or equal to the threshold. We use a conserva¬ 
tive value of 10% of the global peak emission for the threshold, 
which is about ten times higher than the noise level in most of 
the channel maps of all the lines we consider. This conservative 
value is used in order to avoid the association of emission levels 
in one image that would be considered noise in another image. 

Then we multiply each binary channel map image of the two 
line tracers we want to compare, to see if there are regions where 
the two emissions are associated in that particular velocity chan¬ 
nel. The product image would contain pixels with l’s in regions 
where both line tracers have significant emission, and 0’s other¬ 
wise. Thus, adding up all the pixels from the product image, and 
dividing by the total number of pixels, we obtain the fraction of 
the region mapped where the emission from both lines is associ¬ 
ated. An example of this procedure, ap plied to the [C n] 158 pm 
and [C i] 609 pm lines, is shown in Fig. m 

With this method we can estimate the fraction, and where in 
the region mapped, two line tracers are associated at each veloc¬ 
ity channel. The fraction of the region mapped can be compared 
with the correlation coefficient described in Appendix |C| Test 
cases of this are shown in Fig. |D.2| 


Appendix E: [C n] emission associated with the 
three gas phases 

Subtracting the scaled up spectra of [C i] 609 pm, 12 CO J = 
1 —> 0, and the velocity-resolved optical depth r(H i) from the 
[C n] spectra, we obtained a residual [C n] emission that should 
be mostly associated with the ionized hydrogen gas (H n). When 
subtracting the maximum (channel by channel) between the trac¬ 
ers of the H 2 gas (e.g., [C i] 609 pm and 12 CO / = 1 —> 0, or 
[C i] 609 pm and C 18 0 / = 2 —> 1) from the original [C n] spec¬ 
tra, we obtain a second residual [C n] emission that is expected 
to be mostly associated with the neutral (H i) and ionized atomic 
gas (H n). The difference between these residual spectra and the 
original [C n] spectra gives the [C n] emission that is mostly as¬ 
sociated with the H 2 gas. Subtracting the the [C n] emission that 
is mostly associated with the H 2 gas and the first residual spectra 
associated with H n from the original [C n] spectra would lead 
to a third residual [C n] emission that is mostly associated with 
the neutral atomic hydrogen gas, H i. The corresponding column 
densities can be estimated assuming the LTE conditions and Eq. 
[2] The velocity channel maps of the residual [C n] emission as¬ 
sociated to each of the gas ph ases , as estimated with model (2) 
from Table[l] is shown in Fig. |E.l| 




A&A-25020_printer, Online Material p 25 





Fig. C.l. Example of the scatter plots between the pixel values of the velocity-integrated intensity maps of two different line tracers, and the 
corresponding correlation coefficient obtained using Eq. |C.1| 
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Fig. C.2. Example of the scatter plots between the pixel values (Kkms -1 ) of the velocity channel maps of the 12 CO / = 1 —» 0 (X-axis) and 
J = 2-1 (Y-axis) transitions (top left), and the [C I] 370 pm (X-axis) and 13 CO / = 2 —» 1 (Y-axis) lines (top right). The corresponding correlation 
coefficient r at each velocity channel, is shown in the bottom panels. 
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Fig. D.l. Example of the steps in the method to estimate the spatial association between two emission lines. The top panels shows the channel 
maps of the 1 kms -1 integrated intensity (in Kkms -1 ) of the [C II] 158 pm (left) and [C I] 609 pm (right) lines. The contour lines correspond to 
the threshold of 10% of their global peak intensities (136.4 Kkms -1 and 54.6 Kkms -1 for [C n] and [C i], respectively). The middle panels are 
the binary images obtained after applying the intensity threshold to the original channel maps. The bottom panel shows the result of multiplying 
the two binary images, which corresponds to the region where the emission of both [C n] and [C i] lines are associated in a particular velocity 
channel (in this case 15.5-16.5 kms -1 ). 






Fig. D.2. Histograms of the fraction of the region mapped whe re two line tracers are associated at each 1 kms 1 width velocity channel. The 
corresponding correlation coefficient R xy obtained using Eq. |C.l| for each channel map is shown with a continuous line. 
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Fig. E.l. Velocity channel maps at 1 kms 1 width of the [C n] emission (in Kkms 1 ), associated with the H n (top left), HI (top right), and H 2 
( bottom) gas phase, as estimated with model (2) from Table[l] Contours are 20%, 40%, 60%, 80%, and 100% of the respective peak emissions 


Si £ 
























































